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1:  Demonstration  of  Stereo  Mapping  Technology 

1.1  Introduction 

The  focus  of  the  research  over  this  postdoctoral  contract  has  been  in  augmenting  and  refining  the 
specification  of  image  matching  strategies  to  be  employed  in  an  advanced  stereo  mapping  system. 
Initially,  our  research  task  was  to  evaluate  a  preliminary  design  for  a  mapping  system,  and  then 
by  carrying  out  test  implementations  of  elements  of  this  proposed  system,  to  iterate  on  the  design 
process.  As  part  of  a  continuing  mapping  program,  this  phase  of  the  research  would  lead  to  more 
effective  representations  and  improved  matching  strategies  for  surface  reconstruction.  The  test 
implementations  would  build  to  provide  a  foundation  on  which  we  would  develop  an  experimental 
rule-based  stereo  matching  system.  As  part  of  this  experimental  system,  we  undertook  to  develop 
an  interactive  test  facility  to  provide  hand  synthesized  data  for  rule  analysis.  A  change  in  the  scope 
of  the  project  occurred  shortly  after  we  had  begun  the  research,  and  necessitated  a  redirection  in 
the  emphasis;  rather  than  evaluating  or  continuing  the  design,  we  were  to: 

a)  experiment  with  an  existing  matching  system  (showing  its  capabilities 
and  assessing  its  applicability  to  the  planned  system); 

b)  demonstrate  the  design  and  the  utility  of  the  rule-based  approach  to 
surface  inference  from  monocular  information; 

c)  develop  tools  to  support  an  interactive  test  facility. 

1.2  Scope  of  the.  Research 

In  a),  we  undertook  to  apply  an  existing  stereo  mapping  system  [Baker  1981]  to  some  new  imagery. 
This  was  designed  to  demonstrate  its  effectiveness,  to  expose  its  limitations,  and  to  suggest  both 
its  role  in  an  advanced  mapping  system  and  complementary  research  needed  to  improve  its  utility. 
Significant  restructuring  of  the  system  was  called  Tor  in  enabling  it  to  process  this  imagery. 
Details  of  these  changes  are  described  in  section  5,  which  deals  with  the  matching  process.  The 
modifications  have  now  been  implemented,  enabling  the  system  to: 

•  function  on  the  output  of  an  improved  edge  operator  [Marimont  1982]; 

•  use  edge  extent  os  one  of  its  parameters  in  seeking  optimized  correspon¬ 
dence; 

•  exploit  prepared  transform  information  in  processing  images  whose 
epipolar  lines  arc  not  collinear  with  the  scanning  axes  of  the  cameras. 

Separate  funding  sources  supported  research  in  analysis  and  synthesis  of  rules  for  inference  of 
three-dimensional  shape  from  single  images.  This  inference  also  has  application  in  constraining 
search  for  matches  in  stereo  correspondence.  We  will  describe  here  some  results  of  this  research. 
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and,  as  part  of  the  demonstration,  will  show  its  use  in  shape  inference  and  search  constraint  in 
the  domain  of  Orthogonal  Trihedral  Vertices  (OTV’s). 

Our  work  in  developing  tools  has  centered  on: 

a)  a  system  for  the  hand  construction  of  edge  descriptions  from  hard  copy 
imagery; 

b)  an  interactive  system  for  determining  the  transform  to  bring  image 
pairs  into  collinear  epipolar  registration. 

Both  of  these  systems  make  extensive  use  of  interactive  graphics,  and  the  latter  takes  much 
advantage  of  previous  stereo  research  from  our  laboratory  ([Gennery  1980]). 

These  are  the  three  principal  areas  to  be  covered  in  this  report: 

•  experimentation  with  the  automated  mapping  system; 

•  assessment  of  rule  synthesis  from  manual  and  automated  edge  processes; 

•  development  of  test  facility  tools  for  edge  extraction  and  image  registra¬ 
tion. 

1.3  Specific  Tasks 

We  will  describe  results  in  these  areas  of  the  research  through  discussion  of  the  following: 

a)  implementation  of  the  digitising  test  facility; 

b)  the  use  of  image  edge  descriptions  produced  using  this  facility  and 
from  an  automated  process  [Marimont  1982]  in  synthesising  rules 
for  stereo  matching; 

c)  examination  of  Orthogonal  Trihedral  Vertex  (OTV)  inference  rules; 

d)  development  of  a  system  for  registration  of  image  pairs; 

e)  analysis  of  stereo  imagery  with  an  automated  stereo  process. 
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2:  Digitizing  Facility 

2.1  Background 

Our  approach  to  rule  development  begins  with  hand  synthesized  and  some  automatically  generated 
edge  data.  We  have  systems  for  the  automatic  generation  of  edge  data  (i.e.  [Marimont  1982]). 
This  data,  extracted  from  a  sufficently  wide  selection  of  imagery  types,  gives  good  insight  into  the 
current  capabilities  of  automated  processes.  Automated  processes,  however,  are  not  able  presently 
to  give  as  meaningful  a  description  of  an  image  as  we  would  like,  and  have  not  been  designed  to 
provide  the  aggregrated  abstractions  research  systems  ([Lowe  1982])  will  be  soon  supplying.  To 
bridge  this  inadequacy,  we  work  with  both  automatically  generated  data  (the  current  state-of-the- 
art),  and  hand  generated  data  (representative  of  the  next  generation  of  edge  analysis  processes). 
The  hand  generated  data  is  obtained  from  a  manually  operated  digitizing  tablet.  We  have  written 
a  graphics-based  digitizing  and  editing  system  to  run  with  a  GTCO  tablet  in  producing  these 
image  descriptions.  Section  6  lists  instructions  for  the  use  of  this  program  with  the  bitpad  (GTCO 
tablet),  and  section  9  contains  a  source  listing  of  the  program. 

2.2  Imagery  Analysis  with  the  Tablet 

Figure  2-1  below  shows  an  image  pair  of  a  building  complex  (referred  to  as  the  Sacramento 
imagery).  Figure  2-2  shows  the  results  of  bitpad  edge  extraction  on  these  images.  Figure  2-3  shows 
the  results  of  the  Marimont  operator  (Marimont  1982]  on  the  image  pair.  Manually  generated  edge 
data  was  produced  using  this  facility  for  the  analysis  of  rule  synthesis  of  section  3.  It  was  also 
used  to  digitize  the  building  data  of  figures  2-1  for  input  to  the  OTV  inference  process,  as  figure 

3-1. 


Sacramento  Imagery 
Figure  2-1 
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3:  Inference  and  Modelling 

3.1  Introductio.% 

We  have  been  addressing  inference  of  matching  rules  and  the  use  of  model-based  analysis  both 
with  theoretical  analyses  and  with  hand  and  automated  analyses  of  specific  matching  strategies; 
the  latter  applied  to  both  real  and  synthesized  imagery  examples. 

We  have  obtained  extended  edge  data  from  hand  and  automated  processing  for  vr  >  in  synthesis 
of  matching  rules.  Results  from  earlier  work  on  OTV  analysis  (orthogonal  tr  Jr.  I  vertices) 
have  been  exploited,  and  a  demonstration  of  the  use  of  such  a  rule  for  shape  ^rence  and  in 
constraining  search  for  correspondence  has  been  prepared.  Other  funding  has  st  -irted  research 
in  representation  of  generic  structures.  We  have  taken  examples  from  this  moc  **  research  to 
produce  ground  and  aerial  views  of  a  building  complex,  and  have  used  this,  as  w  '  other  data, 
in  rule  synthesis. 

3.2  Geometric  Modelling  and  Computer  Vision 

3.2.1  —  Modelling,  prediction  and  interpretation 

Of  course,  one  of  the  primary  goals  of  research  in  computer  vision  is  the  development  of  systems 
that  can  recognize  and  locate  objects  in  images.  In  order  to  identify  such  an  object,  it  is  clearly 
necessary  to  have  some  description  of  its  characteristics  that  can  be  detected  in  an  image.  A 
representation  of  an  object  is  the  form  this  description  takes. 

One  approach  to  representation  is  to  provide  the  system  with  three-dimensional  models  of  objects. 
Rotation  of  these  models  will  allow  objects  to  be  observed,  conceptually,  from  differing  viewpoints. 
If  parameters  in  a  particular  model  are  allowed  to  vary  it  is  possible  to  have  that  single  model 
represent  a  whole  class  of  objects;  constraining  the  parameters  functions  to  delimit  sub-classes. 
Further  model  manipulations,  such  as  partitioning  and  projection,  can  be  used  to  aid  in  mapping 
model  to  imagery  data.  The  information  contained  in  such  object  models  may  be  used  to  determine 
possible  interpretations  of  image  features  (e.g.  ,  edges,  ribbons,  corners)  and  to  provide  feedback 
to  predict  the  locations  of  such  features  in  an  image. 

ACRONYM  [Brooks  1981]  is  a  three-dimensional  rule-based  modelling/vision  system  developed 
here  at  Stanford  that  provides,  among  other  things,  such  feature  prediction,  model  manipulation, 
and  image  interpretation.  The  rule-base  operates  on  the  models  and  on  the  sensed  data  to 
accomplish  scene  interpretation.  Such  a  rule-based  approach  has  been  shown  to  be  an  effective 
form  for  constraint  and  search  implementation,  and  allows  easy  modification  and  addition  of  new 
rules  without  the  need  of  altering  the  underlying  code. 

Our  group’s  intention  over  the  next  few  years  is  to  build  a  rule-based  stereo  system  operating 
vithin  ACRONYM  whose  functioning  will  include  model-based  prediction.  Working  toward  this, 


6 


Inference  Rules  %S.2.S 


we  have  been  carrying  out  experiments  on  scene  inference  and  model-based  prediction  that  will 
lead  to  a  repertoire  of  stereo  matching  rules. 

3.2.2  -  Models  and  stereo  matching 

One  of  the  major  difficulties  in  determining  stereo  correspondence  is  in  dealing  with  the  large  num¬ 
ber  of  matches  that  are  possible.  Solution  is  generally  found  by  search  through  a  large  parameter 
space,  where  possible  correspondences  are  limited  by  geometric  or  photometric  constraints.  Search 
can  be  reduced  even  more  dramatically  by  endowing  the  matcher  with  broad  domain  specific 
and  domain  independent  knowledge.  Such  knowledge  can  be  rule-based  and  model-based.  Our 
proposition  here  is  that  the  three-dimensional  information  in  object  models,  along  with  inference 
and  prediction  mechanisms,  can  be  used  to  interpret  features  in  image  pairs.  These  interpreta¬ 
tions  can  then  be  used  as  filters  to  constrain  the  matching.  We  demonstrate  this  notion  with  the 
example  of  Orthogonal  Trihedral  Vertices,  often  referred  to  as  cube  corners  or  OTVs.  Other  rules 
synthesized  from  analysis  of  both  manually  extracted  and  automated  edge  processes  follow. 

The  work  on  cube  corners  points  to  additional  usefulness  for  a  model-based  approach.  OTV 
orientation  analysis  (from  matches  across  pairs  of  views)  yields  almost  complete  solution  for  camera 
parrnameters;  constraints  on  sizes  (again,  from  rules  and  models)  could  complete  the  camera 
solution.  Hut  the  orientation  information  yielded  by  a  match  of  a  pair  of  vertices  is  valid  only  if 
the  vertex  is  a  cube  corner;  thus  it  is  necessary  to  be  able  to  distinguish  between  vertices  that  are 
cube  corners  and  those  that  arc  not.  If  the  models  contain  sufficient  information  to  identify  cube 
corners,  then  the  problem  of  determining  cube  corners  independently  of  the  identification  process 
is  eliminated.  In  fact,  both  the  search  for  cube  corners  and  the  search  for  identification  are  likely 
to  be  reduced  when  they  are  combined. 

3.2.3  -  OTV  rule-baaed  analysis 
OTV  theory 

In  cultural  scenes,  we  find  a  large  number  of  interior  and  exterior  corners  of  cubes  typically  when 
two  walls  at  right  angles  meet  the  roof  or  the  floor.  The  importance  of  utilizing  this  common 
structural  element  the  Orthogonal  Trihedral  Vertex  (OTV)  has  been  emphasized  ear!ier[Liebes 
1981)  .Since  they  provide  a  very  tight  constraint  -  the  three  edges  are  mutually  orthogonal  in 
space  it  is  possible  to  calculate  the  three  dimensional  orientations  given  the  projections  in  the 
image.  This  can  be  done  for  both  orthographic  and  perspective  viewing. 

If  the  eye  is  assumed  to  be  focussed  on  the  vertex  of  the  cube  corner,  perspective  can  be  ignored 
and  the  projection  of  a  cube  corner  in  XY Z  space  will  simply  be  its  orthogonal  projection  on  the 
XY  plane.  Suppose  that  some  3  -  star  has  angles  between  its  rays  a,  b  and  r  and  also  that  the 
rays  are  represented  by  the  unit  vectors  w,t  ii2,  t>3.  We  are  interested  in  detecting  whether  there 
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are  three  vectors  in  XY Z  space,  which  are  mutually  orthogonal  and  project,  respectively  to  t>i, 

V2,  V3. 

Since  projection  is  accomplished  by  dropping  the  a-component,  any  3  such  vectors  must  be  of  the 
form  v\  +  Xi.?,  t>2  +  ^2*,  and  v3  +  X3a  where  z  is  the  unit  vector  in  the  ^-direction. 

Requiring  mutual  orthogonality  implies  that  the  dot  products  of  these  vectors  in  pairs  be  zero. 
From  these  conditions  and  some  simple  manipulations  we  can  calculate  the  formulas  for 


K'=±\j- 


(cosa)(cosc) 
(cos  b) 


"f 


(cos  a)(cos  b) 
(cos  c) 


(cos  c)(cos  b) 
(cos  a) 


Hence  solutions  exist  if 

a)  cosa,  cosb,  cosc  are  all  non-zero  and 

b)  either  one  or  three  of  cosa,  cosb,  cosc  are  negative,  so  that  the  quan¬ 
tities  under  the  square  root  sign  are  positive. 

These  results  were  first  derived  in  [Perkins  1968]. 

Thus  we  have  a  way  of  both  eliminating  false  candidates  for  being  OTVs  and  finding  the  3-D 
orientations  of  valid  OTV’s.  This  algorithm  has  been  implemented  and  run  on  data  from  the 
digitizing  tablet. 

OTV  with/ from  models 

Our  analysis  begins  on  both  images,  processing  bottom  up  on  the  two  images  separately.  As  the 
rule  system  identifies  likely  OTVs  in  images  (from  its  models),  it  proceeds  to  match  them.  The 
system  should  already  have  a  tentative  identification  of  the  buildings  containing  the  OTVs,  so 
there  should  be  relatively  few  possible  matches  at  this  point.  Only  OTVs  that  could  be  the  same 
point  on  the  same  object  need  be  compared.  The  analysis  results  in  depths  of  matched  objects, 
for  all  those  objects  having  OTVs. 

This  requires  that  the  modelling  system  handle  point  elements,  and  that  it  include  both: 

•  inferring  OTVs  from  models  (volumes); 

•  accessing  OTVs  stored  explicitly  with  the  models. 
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S.S  Rule  Synthe.no 

3.3.1  -  Inference  n«lcs 

We  continue  with  the  development  of  inference  rules.  This  work  is  a  logical  extension  of  previous 
work  [Binford  1981,  Lowe  1982]  done  at  Stanford  in  developing  mles  for  inferring  surface  infor¬ 
mation  from  a  single  view.  General  assumptions  about  illumination,  object  geometry,  the  imaging 
process  etc.  have  been  used  to  derive  rules  for  making  specific  inferences.  For  Btereo  vision  Arnold 
and  Binford  [Arnold  1980]  have  developed  conditions  on  correspondence  of  edge  and  surface  in¬ 
tervals.  We  divide  our  rules  into  two  categories:  monocular  rules,  which  enable  surface  inference 
from  a  single  view;  and  stereo  rules,  which  facilitate  cross-image  matching. 

3.3.2  —  Monocular  and  stereo  rules. 

1.  Monocular  rules  -  Rules  which  have  been  developed  for  inferences  from 
monocular  views  can  be  utilised  to  provide  a  partial  3-dimensional 
interpretation  which  directs  search  in  the  second  view.  This  category 
includes  the  rule  for  interpretation  of  Orthogonal  Trihedral  Vertices. 

Another  example  is  the  T-junction  rule  [Binford  1981]  which  states  that 
‘In  absence  of  evidence  to  the  contrary,  the  stem  of  a  T  is  not  nearer 
than  the  top,  i.e.  is  coincident  in  space  or  further  away'.  Application 
of  this  rule  gives  a  set  of  nearer/farther  relations.  A  hypothesized 
correspondence  of  edges  which  leads  to  inconsistent  conclusions  from 
the  two  views  can  be  pruned  from  the  search. 

An  image  line  which  is  straight  must  be  the  image  of  a  straight  space 
curve  unless  the  curve  is  planar  and  the  observer  is  coincidentally 
aligned  with  the  plane  of  curvature.  This  enables  us  to  dismiss  cor¬ 
respondences  between  straight  edges  in  one  view  and  curves  in  the  other 
view.  If  two  image  curves  are  projectively  consistent  with  parallel,  we 
assume  they  are  images  of  curves  which  are  parallel  in  space.  That  im¬ 
plies  that  their  images  in  the  other  view  would  be  parallel  i.e.  parallels 
map  to  parallels. 

As  these  examples  illustrate,  most  of  the  rules  in  [Binford  1981,  Lowe 
1982]  and  others  developed  by  Malik  and  Binford  have  as  direct  corol¬ 
laries  stereo  rules  for  checking  the  legality  of  a  match.  They  can  even 
direct  the  search  process. 


] 
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2.  Stereo  Rules  -  these  are  rules  which  have  been  derived  from  the  stereo 
imaging  process,  and  are  a  function  of  the  imaging  geometry. 

An  example  rule  in  this  class,  which  has  long  been  used  for  finding  stereo 
correspondences,  is  the  epipolars  rule  -  corresponding  points  must  lie  on 
corresponding  epipolar  lines.  These  rules  have  inherently  no  monocular 
analogs.  Here  are  a  few: 

a)  Horizontal  planes  in  one  view  get  mapped  to  horizon¬ 
tal  planes  in  the  other  view. 

b)  Use  of  projective  and  quasi- projective  invariants.  This 
has  not  been  examined  in  detail.  Duda  and  Hart(Duda 
1973]  devote  a  chapter  to  this  topic  which  has  not 
really  been  exploited  in  stereo  work. 

c)  Conditions  on  correspondence  of  edges  and  surface 
intervals[Arnold  1980]. 

d)  Surface  Occlusion  rules: 

Surfaces  visible  in  one  view  can  be  occluded  in  the 
other  view.  We  arc  interested  in  the  conditions  when 
this  takes  place.  The  basic  idea  is  that  if  we  cross  a 
surface,  an  obscuration  of'  edge  occurs.  A  left  surface 
visible  in  a  right  view  is  visible  in  the  left  view  un¬ 
less  there  is  obscuration  by  a  tall  object.  Similarly  a 
right  surface  visible  in  a  left  view  will  be  seen  unless 
obscured  by  a  tall  object.  These  surface-obscuration 
rules  can  be  formalized  by  the  cross-product  rule: 


For  the  hypothesized  edge  match  cj  with  f\  and  ej 
with  fs,  we  compute  the  Z-component  of  the  vec¬ 
tor  cross-product  in  the  left  image  pair  and  the  right 
image  pair.  If  the  z-components  have  opposite  signs, 
ve  are  seeing  opposite  sides  of  the  surface.  That  im¬ 
plies  that  the  object  is  not  opaque. 
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3.3.3  -  Use  of  inference  rule*  in  »  test  analytii 

Our  preliminary  results  indicate  good  potential  for  the  success  of  this  approach.  On  hand  simula¬ 
tions  with  line  drawings  of  stereo  pairs,  the  rules  helped  narrow  down  the  choices  considerably. 

Consider  the  imagery  shown  in  figures  3-1  and  3-2  Figure  3-1  is  the  right  view  and  3-2  the  left 
view.  Vertices  1,  2,  3  are  orthogonal  trihedral  vertices.  Using  the  formulae  developed  earlier,  we 
can  find  the  3-D  orientations  of  the  edge  vectors.  These  can  be  matched  with  the  3-D  orientations 
of  l',  2',  3'  to  obtain  a  registering  of  these  vertices  when  combined  with  the  epipolar  constraint. 
All  OTV’s  in  one  view  need  not  be  visible  in  the  other  ue.  4'.  Of  the  monocular  constraints,  the 
other  major  constraints  which  can  be  seen  here  are  the  T-junction  rule  and  the  parallels  rule.  In 
figure  3-2  edge  5  is  behind  edge  6.  Edges  7  and  8  are  parallel  and  so  are  7'  and  8'.  A  match  of  8 
with  9'  would  not  be  accepted.  Surfaces  S\  and  S?  are  both  horizontal  planes  (as  can  be  deduced 
from  the  OTV  analysis)  and  can  be  matched.  Surface  S3  is  not  horizontal.  Here  of  course,  this 
does  not  provide  any  new  information.  Surface  S4  being  a  left  face  in  a  left  view  is  not  guaranteed 
to  be  visible  in  the  other  view  -  as  in  fact  it  is.  The  cross-product  rule  could  be  used  to  di-  miss 
a  match  between  10  and  107. 


Description  of  Left  Image  of  Stereo  Pair 
Figure  3-1 
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4:  Image  Registration 

4-1  Introduction  to  Epipolar  Geometry 


The  search  process  in  automated  stereo  napping  can  be  greatly  restricted,  and  computation 
times  significantly  reduced,  if  information  is  available  relating  the  relative  camera  geometries  of 
a  stereo  pair.  Often  this  information  is  available  in  the  reconnaisance  data  (or  at  least  a  rough 
approximation  to  it).  Other  manual  and  automated  schemes  have  been  devised  to  provide  the 
information  when  it  is  not  present  with  the  imagery  ([Gennery  1980],  [Hallert  I960]).  This  camera 
geometry  information  allows  establishing  epipolar  correspondence  of  lines  across  images.  When 
this  has  been  done,  search  in  one  image  for  match  points  of  a  feature  in  the  other  image  can  be 
constrained  along  a  single  vector.  More  generally,  any  features  lying  along  a  particular  vector  in 
the  one  image  may  be  found  along  a  single  vector  in  the  other  image.  These  image  plane  vectors 
are  termed  epipolar  lines.  Corresponding  vectors  are  termed  corresponding  or  conjugate  epipolar 
lines. 

In  this  section  we  detail  an  algorithm  for  determining  conjugate  epipolar  lines  in  a  set  of  imagery 
for  which  such  camera  geometry  information  is  not  explicitly  available.  Here,  we  rely  upon  an 
operator  to  select  corresponding  points  in  the  two  images.  The  system  automatically  improves 
the  resolution  of  the  correspondence  through  Fourier  interpolation  over  a  match  window  {Gennery 
1980].  The  set  of  such  points  is  taken  by  an  automated  camera  solver  to  produce  the  needed 
geometric  information.  This  point  selection  is  done  with  pan/zoom  cursor  control  on  a  graphics 
device.  If  the  camera  information  is  available  (as,  for  example,  from  reconnaisance  data),  then  the 
point  matching  phase  may  be  omitted  (although  this  provision  has  not  been  enabled  in  the  current 
system).  Equally,  rough  camera  geometry  information,  if  available,  may  be  used  to  partially 
automate  the  point  selection  phase,  although  again  this  is  not  implemented  here.  [Gennery  1980] 
and  [Moravec  1980]  have  implemented  totally  automatic  camera  solvers  in  their  stereo  matching 
systems.  Our  next  improvement  to  this  system  will  be  to  incorporate  the  image  sampling  and 
feature  matching  of  the  [Gennery  1980]  system,  removing  the  need  for  manual  point  selection. 
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4‘2  Glossary  of  Term« 

Given  two  cameras  C\  and  Cg  with  origins  6\  and  0g  and  focal  planes  Pi  and  Pg,  we  call: 

Bundle  of  lines:  A  family  of  lines  containing  a  common  point. 

Epipolar  coordinates  of  a  point:  The  number  of  the  epipolar  line  it  belongs  to,  and  its  distance  to 
a  fixed  reference,  like  the  epipole  if  it  exists. 

Epipolar  direction:  The  direction  of  all  epipolar  lines  if  they  are  parallel. 

Epipolar  line:  The  intersection  of  an  epipolar  plane  with  a  focal  plane.  Alternate  definition:  the 
image  in  one  camera  of  the  pre- image  of  a  point  in  the  other  camera’s  focal  plane. 

Epipolar  plane:  Any  plane  containing  the  two  camera  centers  0\  and  0g. 

Epipolar  space:  A  space  where  the  coordinates  are  the  epipolar  coordinates.  In  this  space,  a 
horizontal  line  is  an  epipolar  line,  and  the  epipole,  if  it  exists,  is  a  whole  vertical 
line. 

Epipole:  The  intersection,  if  it  exists,  of  all  epipolar  lines  in  a  focal  plane. 

Conjugate  epipolar  lines:  the  intersections  of  an  epipolar  plane  with  the  two  focal  planes. 


H 


Epipolar  Geometry 
Figure  4-1 
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4-3  Background  Theory 

Given  two  stereo  images  and  Pt  (the  content  of  the  focal  planes  P\  and  P2  of  the  cameras),  and 
two  lines  L\  and  Lj  contained  in  P\  and  P 2,  in  general  every  point  of  L\  maps  to  a  line  segment  i; 
Pj,  and  there  is  no  particular  relationship  between  the  line  segments  mapping  to  different  point- 

Fundamental  property: 

Every  point  in  L\  will  map  to  a  line  segment  contained  in  the  same  line  if  and  only  if  L\  an- 
L 2  are  a  pair  of  conjugate  epipolar  lines. 

Epipolar  Geometry: 

The  family  of  epipolar  lines  in  a  focal  plane  is: 

Either 

a)  a  set  of  parallel  lines  having  a  common  epipolar  direction, 
or 

b)  a  bundle  of  lines,  the  intersection  of  which  is  the  epipole. 

4-4  Requirements  of  the  System 

Given  a  pair  of  stereo  images,  we  want  to: 

1)  identify  the  kind  of  epipolar  geometry  present  in  the  images; 

2)  explicitly  show  the  epipolar  lines  belonging  to  each  image; 

3)  for  each  image,  compute  the  parameters  which  relate  the  original  coor¬ 
dinates  to  the  epipolar  coordinates; 

4)  construct  the  image  transforms  in  epipolar  space. 

Prior  to  these  4  steps,  we  will  need  to  solve  for  the  cameras,  that  is,  to  determine  the  5  parameters 
describing  their  relative  orientation.  A  procedure  developed  at  Stanford  [Gennery  1980]  is  used 
for  this.  We  proceed  as  above,  using  simple  analytic  geometry  for  our  calculations. 

4-5  Algorithm  Used 

4.5.1  -  Camera  registration  output 

Each  camera  is  viewed  as  a  referential  (0,,  x,  y,  z),  i  €  {1, 2}.  The  registration  procedure  yields 
azimuth,  elevation,  pan,  tilt,  and  roll  of  one  camera  with  respect  to  the  other.  From  there,  ”/e 
compute: 


Eptpolar  Geometry  §^-5.i 


Too  view  o/a  situation  with  only  one 

epipoie 

Figure  4-2 

IM 
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\ 

Too  view  of  a  situation  with  epipoie j 
Figure  4-3 
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1)  the  rotation  matrix  R  between  the  two  referentials: 


2)  the  translation  unit  vector  t,  the  components  of  which  are: 


0  |in  the  base  0\xyz, 


H  Jin  the  base  0? xyz 


Note  that  the  magnitude  of  the  translation  vector  cannot  be  determined  from  a  pair  of  images. 
4.5.2  -  Epipolar  geometry  determination 

The  focal  planes  Pt  are  planes  parallel  to  0,xy,  intersecting  0,z  at  z  —  fit  the  focal  distance. 
There  is  an  epipole  in  plane  P,  if  and  only  if  the  translation  vector  intersects  this  plane,  that  is, 
if  and  only  if  its  third  component  is  not  sero. 

We  thus  determine  the  case  we  are  working  with: 


if  7  5^  0  and  v  ^  0,  there 
if  7  0  and  v  =  0,  there 

if  7  =  0  and  u  ^  0,  there 
if  7  =  0  and  u  =  0,  there 


are  two  epipoles: 

CASE  1 

is  one  epipole  E\: 

CASE  2 

is  one  epipole  E 2: 

CASE  3 

are  no  epipoles: 

CASE  4 

Remarks; 

a)  7  =  0  is  replaced  in  the  code  by  |7|  <  threshold,  where  threshold  is 
chosen  as  a  function  of  the  arithmetic  precision  of  the  machine:  if  we 
had  infinite  precision,  then  we  could  consider  every  case  as  being  case 
1.  Here  threshold  =  .0001  was  found  to  be  a  good  estimate. 

b)  Most  image  pairs  will  belong  to  the  first  case,  with  7  and  v  of  the  order 
of  .1.  The  epipoles  exist,  are  outside  the  picture  frame,  and,  for  the 
images  worked  with  to  date,  tend  to  be  at  a  distance  of  about  10  times 
the  picture  dimension. 
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4.5.3  —  Epipolea  and  epipolar  directions 

If  the  epipole  Ei  exists,  it  is  the  extremity  of  the  vector  collinear  to  the  translation  vector,  with  a 
third  component  equal  to  /,-.  If  it  does  not,  then  the  translation  vector  is  the  epipolar  direction. 
Hence: 


Case  1: 

Case  2: 

Case  3: 

Case  4: 

4.5.4  —  Epipolar  line  calculation 

4.5.5  -  CASE  Is  two  epipoles 
a)  theory 

Let  M\(x\,y\,Z\)  be  a  point  in  P\.  EXMX  defines  an  epipolar  line  in  Px,  and  the  corresponding 
E2M3  defines  the  conjugate  epipolar  line  in  P2.  The  plane  (f?i,0j,<?2>  £2,  A/,,  M2)  contains  the 
translation  vector  t  and  EXM\.  Its  normal  is  E\M\  X  t-  The  normal  of  P-i  is  02*.  Hence  the 
intersection  of  the  two  planes  is  given  by  the  vector: 


O^z  X  (£1  M\  X  t)  =  (02a  •  t)Ei  M\  —  (022  ■  ExMx)t 

and  £2A/2  is  collinear  to  this  vector.  Suppose  that  in  9xxyz,  EXMX :  (zt,yi,0).  In  terms  of 
components  in  02zyz: 

0f\\  fx\\  /r,,ii  +  riayA 

,  t  =  EiMt  =  I  y\  J  =  (  r2,zi  +  r22y,  1 

\vj  \z\)  +  *-32Vl/ 


E2  M2  is  collinear  to 


'vx\  -  Xz^X  _  / xi(uru  -  Xr3I) +  y x(vria  -  Xr32)X 
^y'l  -  HA)  -  /W31)  +  V\(vrn  -  fir 32) J 
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If  we  let  A  be  the  matrix: 


(vrn-  XrSi 

V^ai  —  Mr3i 


ur12  -  \ria 
vru  -  fir 3j 


Then  we  can  write  EaMa  =  AE\ Mi  where  E\M\  is  in  base  Ojxy  and  E2M2  is  in  base  02xy. 


P. 


Pi- 


Case  1 
Figure  4-4 

b)  algo rithm 

Let  N i  be  the  number  of  epipolar  lines  that  we  want  to  determine.  Each  epipolar  line  is  uniquely 
determined  by  the  angle  it  makes  with  the  z-axis.  Let  0  be  this  angle.  Given  k,  the  epipolar  line 
number,  0  <  k  <  nl  —  1,  how  can  we  determine  01  If  0 o  and  0\  are  the  lower  and  upper  limits 
between  which  0  is  allowed  to  vary,  and  02  =  — ^ ,  then  we  will  choose  the  middle  of  each 
interval:  0  —  0$  +  (k  +  .5)02  But  what  are  Oq  and  0i?  We  have  to  distinguish  between  three  cases: 

•  The  epipole  is  in  the  picture  (very  unlikely).  Then  0  can  vary  between 
0  and  2x  radians; 
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•  The  epipole  is  outside  the  image:  there  is  a  minimum  and  maximum 
angle  under  which  the  image  is  seen  from  this  point.  If  we  choose  these 
angles  in  [0,2irj,  then  most  'f  the  time  every  $  G  [#o>0i]  will  define  a 
valid  epipolar  line  in  Pi; 

•  The  exception  from  above  is  when  the  epipole  is  left  of  the  image  but  on 
a  same  vertical  level:  then  [tfo.^il  cannot  be  connected  and  still  included 
in  [0,2*].  In  this  case  we  will  choose  the  angles  in  [^,  5]. 


Then  L\  will  be  defined  by  the  point  E\  and  the  vector  V](co80,  stn0),  and  L2  is  defined  by  E2 
and  V2  =  AV\. 


4.5.6  -  CASE  2:  one  epipole  E\ 


a)  theory 


Given  an  epipolar  line  ,  we  already  know  that  the  corresponding  epipolar  line  L2  in  P2  is 

coliinear  to  the  vector  t  =  V2(\,fi).  Hence  we  just  need  to  find  a  point  belonging  to  L2.  Clearly, 
ha  is  the  intersection  of  P2  with  the  plane  (E 1,  M\,0 2).  Thus  any  line  contained  in  this  plane  will 
intersect  P2  at  a  point  contained  in  L2.  In  particular,  consider  the  parallel  to  L]  driven  through 
0j.  It  intersects  P2,  thus  L2,  at  M3  such  that,  in  base  02xyz: 


TiiXi  +  r22yi  1, 
ra\X\  +  r3  3yiJ 


hence  62M3  = 


/  r„  ' 

/  •'^'31*1 +MaVi) 

a,  lL».t,->+.r»»K> 


fl 


b)  algorithm 

In  the  same  way  as  in  case  1,  we  define  L\(E\,  V\ )  where  V\(x\  =  co»0 \,x2  =  sin&2).  Then  La  is 
defined  by  {M3,  V2),  where  the  coordinates  of  A/j  are 


r  11X1  +  r  1 2  7/ 1 
ra  i*\  +  r3iVi 


h 


r2ixi  4  r;aV! 
r3J  xl  +  r32'Jl 


1 

I 

I 
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Case  2 
Figure  4-5 


Let  an  epipolar  line  in  Py  be  defined  by  the  translation  vector  t  =  Vy  and  by  a  point  My  we 
pick.  In  P2,  Li  goes  through  £2  and  is  collinear  to  a  vector  Hi  Mi,  intersection  of  Pi  with  the 
plane  {62,  E2,  My).  This  plane  is  orthogonal  to  02My  X  t  and  P2  is  orthogonal  to  02z.  Hence  the 
intersection  is  collinear  to  02z  X  [02M\  X  t) 

since  62 My  =  0201  +  0y  My 
=  kt  4-  OyMy , 

02My  X  t  =  OyMy  X  t, 

and  $2 Z  X  [Q2My  X  t)  =  (02Z  •  l)0yMy  —  {32Z  ■  OyMy)t. 
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Suppose  0i  Mi  :  (xi,yt,/i)  in  B\xyz.  Then  in  O^xyz: 


E2M2  is  thus  collinear  to: 


(full  +  r12y,  +  ri3/j 
f2iZi  +  r2  2V1  +  *"33/1 
^31 2=1  +r32yi  +  r33/l 


fvx\  -  Xy,\  _  /  xt (i/ru  -  Xr31) +  yi(t/r,2  -  Xr32)  +  /i(rj3  -  Xr33)  \ 
\uy\  -  pz'i)  \ii(t'r2i  -  pr3i)  +  yi(t /r22  -  /ir32)  +  fi(ur23  -  ^33)) 


Hence,  if  we  let  A  be  the  same  matrix  as  in  CASE  1  and  OF3  be  the  offset  matrix: 


/Vr,3  -  Xr33\ 

\i/r 23  ~  Mr33/ 


Then  we  have:  E2M3  =  AOtMi  +  OF 3 

Where  OtMi  is  in  base  0\xy  and  E2M2  is  in  base  02xy 

b)  algorithm 

Now,  how  do  we  pick  M\  in  the  first  place?  We  want  a  set  of  N\  equally  spaced  epipolar  lines, 
and  it  appears  convenient  to  pick  points  on  the  axes.  If  the  epipolar  lines  are  more  horizontal, 
or  the  image  stretched  in  height,  then  we  will  pick  N\  equally  spaced  points  on  the  vertical  axis, 
suitably  located  to  cover  the  entire  image.  If  the  epipolar  lines  are  more  vertical  or  the  image 
more  stretched  in  width,  then  we  pick  them  on  the  horizontal  axis.  Let  LZ,LV  be  the  picture 


dimensions  and 

(V 

*,Vy) 

the  epipolar  direction: 

If  VxLy 

< 

VyLZ, 

and 

v. 

<  0, 

we 

pick 

Vi 

=r 

+  ^v)( 

fC+0.5\ 

If  VtLy 

< 

VyLX, 

and 

v* 

>  0, 

we 

pick 

y\  = 

=  (LZ\ 

Iv^l  +  ^v)( 

/C+0.5A 

t- 

If  VvLz 

< 

VZLy, 

and 

Vy 

<  0, 

we 

pick 

= 

+  £*)( 

rc+0.5  ^ 

If  VyLx 

< 

V*Ly, 

and 

V'v 

>  0, 

wc 

pick 

X!  = 

=  (Ly 

l%l  +  M(*£rL)- 

Then  wc  proceed  as  indicated  abo/e: 

=  t)  is  matched  with  L^(E2,  V2  =  F2M2). 
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4.5.8  -  CASE  4:  no  epipoles 


a)  theory 

Let  an  epipolar  line  in  P\  be  defined  by  the  translation  vector  t  =  V\  and  by  a  point  Mj  we  pick. 
In  P^,  L 2  goes  through  the  image  of  M\,  that  is  the  extremity  of  a  vector  collinear  to  0\Mi  and 
whose  third  component  is  S2 ■  In  base  0xxyi,6i  M\  :  In  base  02XV*: 


OjMi 


(r  n*t  +ri2Vi  +r13/, 
+  r22Vl  +  rj3/i 
r31xl  +  ’’321/1  +  »,33/l 


rll  zl  +rnVi+rn/i 
^31*1  +^32  Vl  +  r33jT» 
*31  *  I  +  *33  Vl  +  *33/1 
*3i  *  1  +*33yi+*33/i 

ft 


) 


b)  algorithm 


I 


We  pick  points  E 1  in  the  same  manner  as  in  case  3.  Then  we  calculate  O2M2  as  indicated  above 
and  we  match  L\{E\,t )  with  L2(E2,t). 
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4-6  Transform  Parameters 

We  want  to  output  the  minimum  necessary  parameters  to  do  the  following  transformation: 
given  the  coordinates  of  one  point  in  one  image,  find  its  coordinates  in  epipolar  space. 

We  treat  the  problem  only  in  CASE  1,  which  is  the  most  common  case  encountered.  For  this  case, 
we  will  need: 

•  the  coordinates  of  the  epipoles; 

•  the  number  of  epipolar  lines; 

•  the  minimum  and  maximum  angle  0q  and  0\  under  which  picture  1  is 
viewed; 

•  the  minimum  distances  D\min  and  Djmt'n  of  the  epipoles  to  the 
images; 

•  the  matrix  R,  or  rather  its  inverse  R-1. 
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For  example,  the  images  of  figure  5-6  have  transform  parameters  as  follows: 


Case  =  1 


N  =  450 
En  =  356755.3 
5,3  =  166509.8 
Eai  =  -7765.348 
E22  =  -3509.268 
0O  =  3.578705 
0,  =  3.577792 
D\min  =  393249.0 
Damin  —  8521.478 
RTS  =  430.7463 
RTl  =  2375.517 
R2i  =  215.9756 
R22  =  1075.811 


4-7  Example  of  Epipolar  Registration  and  Tranoformation 


Figure  4-6  shows  a  stereo  pair  of  a  building  complex.  Figure  4-7  has  this  pair  superpositioned  with 
a  net  of  corresponding  epipolar  lines.  Figure  4-8  shows  the  imagery  transformed  such  that  epipolar 
lines  arc  horizontal  in  the  image,  and  conjugate  epipolar  lines  have  the  same  row  coordinate. 
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5:  Analysis  of  Automated  Stereo  Mapping 

5.1  Background 

Results  from  our  laboratory  over  the  past  few  years  [Quam  1971,  Hannah  1974,  Moravec  1980, 
Gennery  1980,  Arnold  1980,  Baker  1981,  Arnold  1983],  have  demonstrated  the  possibilities  of  both 
area- based  and  feature- based  stereo  matching. 

Area-based  stereo  matching  uses  windowing  mechanisms  to  isolate  parts  of  two  images  for  cross¬ 
correlation.  Feature-based  stereo  matching  uses  two-dimensional  convolution  operators  (and  per¬ 
haps  grouping  operators)  to  reduce  an  image  to  a  depiction  of  its  intensity  boundaries,  which 
can  then  be  put  into  correspondence.  Area- based  cross-correlation  techniques  require  distinctive 
texture  within  the  area  of  correlation  for  successful  operation.  They  break  track: 

•  where  there  are  ambiguous  textures  or  featureless  areas  (roofs,  sand 
and  concrete); 

•  where  the  correlation  area  crosses  surface  discontinuities  (at  occlusions 
such  as  buildings,  or  thin  objects  (poles)); 

•  where  depth  is  ill-defined  (such  as  through  trees). 

In  general,  these  systems  break  track  where  there  is  no  local  correlation  (zero  signal,  or  where  two 
images  do  not  correspond)  or  where  the  correlation  is  ambiguous  (where  the  signal  is  repetitive). 
The  systems  must  be  started  manually  and  corrected  when  they  break  track. 

Demands  of  mapping  in  cultural  sites  and  in  locales  with  surface  discontinuity  and  ambiguous 
or  non-existent  texture  make  it  essential  that,  if  area-based  analysis  is  to  be  done,  it  be  done 
in  conjunction  with  feature-based  analysis.  Feature- based  analysis  provides  a  solution  to  many 
of  the  problems  of  correlation.  Principal  among  its  advantages  is  that  it  operates  on  the  most 
discriminablc  parts  of  an  image:  places  that  are  distinctive  in  their  intensity  variation,  and  where 
localization  is  greatest.  These  are  typically  the  boundaries  between  objects  or  between  details  on 
objects,  or  between  objects  and  their  backgrounds.  The  important  point  is  that  the  features  being 
put  into  correspondence  for  depth  estimates  are  the  boundaries  of  objects:  area-based  analysis  is 
at  its  worst  at  object  boundaries,  yet  determining  boundaries  can  be  said  to  be  the  most  important 
part  of  mapping  in  3-space. 


28 


Automated  Mapping  §5.£ 


Some  other  advantages  of  feature- based  analysis  are: 

a)  reduced  combinatorics  —  there  are  fewer  features  than  pixels, 

b)  greater  accuracy  —  features  can,  in  general,  be  positioned  to  sub- pixel 
precision,  while  area  positioning  precision  is  inversely  proportional  to 
window  size,  and  considerably  poorer,  and 

c)  more  realistic  invariance  assumptions  —  area-based  analysis  presup¬ 
poses  that  the  photometric  properties  of  a  scene  are  invariant  to  viewing 
position,  while  feature-based  analysis  works  with  the  assumption  that 
it  is  the  geometric  properties  that  are  invariant  to  viewing  position. 


The  Baker  system  [Baker  1981]  is  the  only  current  system  that  mixes  these  two  matching 
modalities.  We  undertook  in  this  postdoctoral  work  to  demonstrate  the  capabilities  of  this  Bystem 
applied  to  new  cultural  imagery. 


5.2  Baker  System:  before 

In  the  Baker  system,  image  edges  are  the  features  considered  in  stereo  matching.  These  edges 
are  located  at  positions  in  the  image  where  a  change  in  sign  of  second  difference  in  intensity 
occurs.  In  the  initial  case  a  particular  operator  (one  that  was  1  by  7  pixels  in  size)  measured  the 
directional  first  difference  in  intensity  at  each  pixel.  Second  differences  were  computed  from  these, 
and  changes  in  sign  of  these  second  differences  were  used  to  interpolate  zero  crossings  (».e.  peaks 
in  first  difference).  Certain  local  properties  other  than  position  were  measured  and  associated  with 
each  edge  —  contrast,  orientation,  and  intensity  to  either  side  —  and  links  were  kept  to  nearest 
neighbors  above,  below,  and  to  the  sides.  It  is  these  properties  that  define  an  edge  and  provide 
the  basis  for  the  correspondence  process. 

The  correspondence  is  a  search  for  edge  matches  across  images.  With  no  prior  knowledge  of  the 
viewing  situation,  one  could  have  any  edge  in  one  image  matching  any  edge  in  the  other.  The 
combinatorics  of  a  naive  matching  strategy  clearly  could  be  enormous.  A  lot  of  the  analysis  of 
a  matching  strategy  goes  into  bounding  this  combinatorics  and  constraining  the  search  for  edge 
correspondences. 

One  of  the  principal  constraints  used  in  stereo  analysis  may  be  determined  from  a  knowledge  of 
the  relative  attitudes  of  the  cameras.  If  two  rqualivalent  cameras  are  arranged  with  axes  parallel, 
as  shown  in  Figure  5-1,  then  they  can  be  conceived  of  as  sharing  a  single  common  image  plane. 
Any  point  in  the  scene  will  project  to  two  points  on  that  joint  image  plane  (one  through  each  of 
the  two  lens  centers),  the  connection  of  which  will  prrdrce  a  line  parallel  to  the  baseline  between 
the  cameras.  Thus  corresponding  edges  in  the  tv/n  images  muct  lie  along  the  same  line  in  the 
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joint  image  plane.  This  line  is  termed  an  epipolar  line  (see  [Hallert  1960],  and  section  4).  If  the 
baseline  between  the  two  cameras  happens  to  be  parallel  to  the  scanning  axis  of  the  cameras,  then 
the  correspondence  only  need  consider  edges  lying  along  matched  lines  parallel  to  that  axis  in  the 
two  images.  These  lines  are  termed  conjugate.  Figure  5-1  indicates  this  camera  geometry  a 
geometry  which  produces  eollinear  epipolar  lines.  The  algorithm  described  here  assumed  the  stereo 
pair  to  have  such  a  geometry,  and  if  this  was  not  the  case  then  the  appropriate  transformation 
of  the  two  images  had  to  be  made  before  any  processing  could  be  done.  A  less  restrictive  solution 
would  be  to  have  the  correspondence  process  informed  of  the  camera  geometries,  and  have  it  solve 
for  the  more  general  epipolar  situation  as  shown  in  figure  5-2.  Incorporation  of  this  capability  is 
one  of  the  enhancements  we  have  made  to  the  system  over  the  past  few  months.  Section  4  details 
the  analysis  for  this  solution. 


Specialized  Epipolar  Geometry 
Figure  5-1 
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Ep / po / an  linos 


Genera)  Epipolar  geometry 
Figure  5-2 


The  advantage  of  processing  along  epipolar  lines  is  that  the  search  for  correspondence  may  be 
constrained  to  one  dimension,  along  a  single  line-pair.  It  is  then  feasible  to  apply  an  efficient 
optimization  process  to  the  matching  problem.  The  Viterbi  algorithm  [Forney  1973]  is  the 
optimization  process  used  in  this  system.  It  is  a  dynamic  programming  technique  used  extensively 
in  speech  processing,  and  first  used  in  vision  research  in  some  recent  work  at  Control  Data 
Corporation  [Henderson  1979).  An  earlier  use  of  a  dynamic  programming  technique  for  stereo 
matching  is  documented  in  [GimelTarb  1972).  [Baker  1981]  describes  in  detail  the  use  of  the 
Viterbi  algorithm  for  stereo  correspondence. 

Other  geometric  and  photometric  constraints  may  be  introduced  both  in  limiting  search  and 
in  bounding  parameters  for  the  evaluation  function.  [Arnold  1980]  and  [Baker  1981]  provide 
discussions  of  these  constraints  in  determining  edge  correspondences. 
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5.3  Baker  System:  enhancement* 

As  stated,  we  undertook  in  this  postdoctoral  work  to  demonstrate  the  capabilities  of  the  Baker 
system  (Baker  1981)  applied  to  some  cultural  scenes.  Before  carrying  out  these  analyses  we: 

a)  enhanced  the  system  with  a  capability  to  work  with  a  better  edge 
operator  (Marimont  1982); 

b)  enabled  it  to  process  images  that  are  not  graced  with  collinear  epipolar 
geometry  («.e.  most  images); 

c)  introduced  an  additional  correspondence  measure  -  edge  extent. 

To  implement  these  enhancements  required  substantial  redesign  of  the  system,  and  redesign  cycles 
with  the  Marimont  process.  Chosing  useable  data  also  presented  difficulties,  as  the  only  imagery 
available  was  not  of  the  correct  geometry  (see  below).  The  two  image  pairs  initially  chosen  (the 
Sacramento  apartment  complex  and  a  section  of  some  imagery  of  Moffett  Field)  proved,  on  closer 
examination,  to  require  quite  complex  transformation,  and  could  not  be  easily  adjusted  for  epipolar 
processing.  Obtaining  useable  data  then  tied  the  progress  of  this  part  of  the  demonstration  to  the 
progress  of  developing  an  image  registration  system. 

In  general,  to  bring  imagery  data  into  a  properly  transformed  state  could  proceed  in  one  of  two 
ways: 

•  one  could  determine  the  transforms  and  then  modify  the  imagery, 
producing  an  image  pair  having  collinear  epipolar  geometry; 

or 

•  one  could  determine  the  transforms,  and  modify  the  output  of  an  edge 
operator  process  that  functions  over  the  original  imagery. 

The  latter  is  by  far  the  superior  approach,  as  it  avoids  resampling  the  image.  This  approach 
necessitates  incorporating  the  transform  computation  into  the  stereo  system,  to  follow  edge  finding 
and  precede  edge  matching. 

The  second  part  of  the  stereo  system’s  analysis  is  an  intensity  correlation  process.  This  operates 
along  epipolar  lines  as  well,  and  clearly  requires  intensity  information  to  be  accessable  along 
epipolar  lines.  One  solution  to  this  would  be  to  take  the  original  image  pair  and  have  the  correlator 
rotate  and  change  shape,  size,  and  orientation  as  it  moves  around  the  image;  this  is  an  awkward  and 
probably  unnecessary  complication.  An  alternative  would  be  to  access  the  transformed  images, 
sampled  as  accurately  as  possible,  and  do  the  correlation  in  the  rectangular  space  defined  by 
collinear  epipolar  lines.  The  argument  from  edge  accuracy  indicated  that  transforming  edges 
rather  then  resampling  the  image  was  the  way  to  go;  this  argument  from  intensity  correlation 
suggests  that  the  resampled  image  can  to  us.f  il. 
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Another  implementation  detail  supported  this  use  of  both  transformed  edges  and  transformed 
imagery:  it  was  found  that  the  intensity  information  available  from  the  Marimont  process  had 
too  small  a  basis  for  useful  correlation,  and  in  fact,  for  transformed  edges,  had  little  relevance 
for  the  matching  (it  being  measured  not  along  epipolar  lines).  The  transformed  image  had  to 
be  referenced  again  by  the  system  to  obtain  more  significant  intensity  estimates  oriented  along 
epipolar  lines,  and  working  with  the  image  in  epipolar  space  facilitated  this. 

The  philosophy  of  the  stereo  matching  process  here  had  been  to  use  edge  analysis  for,  among 
other  things,  its  higher  accuracy,  and  to  use  intensity  analysis  for  the  continuity  it  provides.  To 
be  consistent  with  this,  we  wanted  to  have  the  highest  possible  accuracy  for  edges  in  epipolar 
space,  and  if  sacrifice  be  needed  for  simplicity,  to  do  it  where  it  least  degraded  the  analysis  - 
in  the  intensity  correlation.  It  is  clear  that  transformed  edges  give  higher  accuracy  than  edges 
from  transformed  images  (detectability  might  not  change  much,  but  localization  is  significantly 
reduced);  and  important  simplifications  could  be  obtained  for  little  loss  by  doing  the  intensity 
correlation  over  the  resampled  image  pair.  This  meant  changes  in  our  plans  for  the  registration 
system:  it  had  to  produce  not  just  transform  information,  but  transformed  images  as  well.  Both 
forms  are  made  available  as  output  from  the  registration  program  described  in  section  4,  and 
the  enhanced  Baker  system  uses  them  both.  (However,  the  original  edge  finder  of  the  stereo 
system  has  not  been  modified  to  use  the  transform  information,  and  must  work  (at  present)  on 
the  transformed  images  with  the  limitations  that  brings.) 

5.4  The  Marimont  Edge  Operator 

The  Marimont  edge  operator  has  greater  detection  and  reliability  than  the  original  Baker  edge 
operator,  and  similar  localization;  earlier  examples  of  its  processing  convinced  us  that  its  output 
would  improve  the  quality  of  our  stereo  reconstruction.  Its  ability  to  track  along  zero  signal 
areas  in  following  zero-crossing  edges  leads  to  more  coherent  image  descriptions.  (Marimont  1982] 
provides  details  of  the  operator’s  functioning.  Roughly,  it  works  by  convolving  an  m  X  to  lateral 
inhibition  function  ofnXn  central  window  with  an  image.  Zero  crossings  in  this  resultant  image 
then  indicate  edges,  and  the  edge  position  is  determined  by  interpolating  over  the  lateral  inhibition 
surface.  Section  7  details  the  content  and  format  of  the  edge  files  produced  by  this  process. 

A  few  unanticipated  problems  became  apparent  once  work  with  the  edges  was  begun.  One  point, 
noted  above,  was  that  the  intensity  information  stored  at  an  edge  (its  left  and  right  boundary 
values)  had  quite  small  support  (a  single  pixel).  This  is  in  contrast  with  the  original  operator  which 
interpolated  for  these  values  in  an  area  3  pixels  wide  and  removed  one  pixel  from  the  determined 
edge  position.  Another  problem  was  that  the  edge  connectivity  produced  by  the  Marimont  system 
can  be  misleading,  as  figure  5-3  shows.  Intensity  significance  was  improved  by  sampling  along 
epipolar  fines  in  the  transformed  images.  The  connectivity  problem  has  not  been  looked  at  yet. 
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Good  connectivity  is  inherently  difficult  to  achieve  with  zero  crossing  operators.  Refinements  to 
the  process  are  being  considered. 


image 
Figure  5-3 


The  introduction  of  edge  extent  as  a  parameter  in  the  dynamic  programming  solution  was  an 
obvious  fallout  from  using  the  Marimont  edges.  Edges  are  output  by  that  process  as  strings,  with 
'■.’-connectedness.  The  maximum  and  minimum  of  some  string,  in  transform  space,  is  a  measure 
of  its  (epipolar)  extent.  Prior  to  the  use  of  this  information  the  only  way  that  global  continuity 
entered  the  analysis  was  through  a  consistency  enforcement  relaxation  process  which  ensured  that 
edges  connected  in  one  view  were  interpreted  as  continuous  in  3-space;  all  matching  measures 
were  i|uitc  local.  W  it  li  t  he  modified  approach,  the  correspondence  measure  is  a  function  of  (among 
other,  more  statistically  based  parameters)  the  ratio  of  edge  extents.  In  particular,  the  likelihood 
of  edge  element  a  in  the  left,  image  matching  edge  element  6  in  the  right  image  depends  on  the 
product  of  the  ratios  of  the  two  upper  extents  (up  from  the  edge  elements)  and  the  two  lower 
extents  (down  from  the  two  odgo  elements). 

When  image  testing  began  with  all  of  the  above  accomplished,  another  problem  became  apparent: 
the  stereo  system,  bound  into  a  machine  architecture  with  a  maximum  of  256K  words  of  memory, 
and  always  tightly  wedged  anyway,  had  grown  with  these  changes  to  the  point  that  only  small 
portions  of  images  could  be  worked  on  at  once.  Thus  came  to  exist  a  windowing  mechanism  within 
the  edge  finding/loading  and  stereo  matching  processes. 


i 
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5.5  Teat  Imagery 

Our  testing  has  been  progressing  on  several  sets  of  imagery:  a  synthetic  image  pair  from  Control 
Data  Corporation,  an  aerial  scene  from  the  Engineering  Topographic  Laboratory,  and  a  building 
scene  of  Sacramento.  The  latter  imagery  is  the  most  relevant  to  the  current  postdoctoral  contract, 
and  we  will  be  demonstrating  the  mapping  system  on  this  data.  This  imagery  may  be  seen  in 
figures  4-6  and  4-8. 
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6:  Bitpad  Instructions 

6.1  Instructions  for  the  Digitizing  Tablet 

6.1.1  -  Introduction 

The  bitpad  program  is  written  in  PASCAL  to  be  run  on  a  DEC  KL-10  processor  under  WAITS 
at  the  Computer  Science  Department  of  Stanford  University.  It  has  an  assembly  code  driver  to 
handle  input  from  the  tablet  over  a  normal  RS-232  port  (listing  included  also),  and  requires  a 
graphics  package  available  locally.  Replacing  this  package  for  an  intended  use  away  from  SAIL 
WAITS  should  be  quite  straightforward. 

Locally,  the  bitpad  may  be  run  by  typing  DO  TAB[S,HI1B].  Figure  6  shows  the  digitizing  tablet 
in  use. 

The  initial  preamble  of  questions  deals  with  graphics  options  and  line  characteristics.  You  will 
usually  want  some  sort  of  graphics  output,  so  answer  yes  to  the  first  question  (“Want  graphics 
output  (y / n )?:” ) .  Graphics  choices  are  DD  or  GOD  (for  later  display).  Say  y  to  either  or  both  of  the 
questions  about  this.  DD  graphics  go  the  Data  Disc  display  while  for  viewing  while  you  construct 
an  edge  file,  while  GOD  graphics  are  for  later  Graphics  Output  Device  device  independent  display. 
If  you  want  little  boxes  drawn  around  the  endpoints  oT  lines,  then  asnwer  y  to  (“Want  boxes  draw 
around  vertices  (y/n)?:”).  Box  size  refers  to  the  size  of  the  box  to  be  drawn  around  vertices.  Small 
is  better  for  entry,  large  makes  it  easier  to  distinguish  vertices  from  line  nodes  if  you  are  trying 
to  do  deletions. 

On  the  first  digitizing  pass  over  an  image  you  will  start  directly  from  the  bitpad;  later  sessions  will 
involve  adding  to  or  editing  an  existing  edge  file.  To  enter  an  old  file  (or  several  old  files),  say  Y 
to  the  question  (“Want  to  read  in  a  file  (y/n)?:”).  This  reads  it  in  and  displays  it.  For  the  moment 
there’s  a  problem  with  reading  strings  into  PASCAL,  so  it  tries  to  read  NEW. VRT[S,HHB].  If 
this  file  doesn’t  exist,  it  will  give  you  the  chance  to  enter  another  filename.,  don’t  panic,  just  type 
in  the  new  name. 

When  that  file  has  been  read  in  and  displayed,  it  lets  you  choose  another  (ad  infinitum).  End  with 
an  answer  of  N  (for  NO)  to  the  above  question. 

6.1.2  -  Setting  screen  coordir.at.s 

Whether  you’ve  read  in  a  file,  or  not,  you  now  set  up  the  bitpad  to  enter  data.  It  asks  for  the 
TTYline  of  the  tablet  (50  is  the  current  line  number  of  the  bitpad  plug).  You  then  define  the 
screen  area  on  the  bitpad:  select  the  four  corners  with  the  crosshairs  (lower  left,  upper  left,  upper 
right,  lower  right)  with  any  button,  and  confirm  your  selection  by  pushing  button  ‘2’  (‘8’  says  no, 
try  for  four  more  points). 


36 


Bitpad  § 6.1.4 


6.1.3  -  Adding  segments  and  vertical 

A  ‘segment’  is  made  up  of  a  start  vertex,  any  number  (even  0)  of  intermediate  points,  and  an  end 
vertex.  Vertices  may  belong  to  any  number  of  ‘segments’.  Adding  vertices/points/segments  is  as 
follows: 

Button:  ‘l’  defines  a  new  vertex  at  the  current  crosshairs  position,  either  starting  a  segment  with 
it,  or  ending  one  there.  ‘2’  locates  the  nearest  vertex  (created  by  ‘1’),  and  begins  or  ends  a  segment 
there.  ‘4’  says  add  this  point  to  the  current  segment  list.  ‘8’  is  a  break.  It  lifts  the  pen.  If  it  is 
followed  by  a  ‘1’  or  ‘2’,  then  that  is  all  it  does.  If  followed  by  a  ‘4’,  it  puts  you  in  the  delete  cycle. 
If  followed  by  another  ‘8’,  you  exit  everything. 

6.1.4  -  Editing  segments 

You  may  enter  the  editor  to  delete  segments  at  any  point.  You  do  this  by  keying  ‘8’  then  ‘4’  in 
the  above  entry  loop. 

To  delete  segments,  first  choose  a  vertex  attached  to  the  segment  with  ‘1’  button.  Then  select 
either  the  other  vertex  (if  there  is  no  ambiguity)  with  a  ‘1’,  or  a  point  somewhere  on  the  segment 
with  a  ‘4’.  You  then  see  it  erased  on  the  screen,  and  you  push  ‘2’  to  confirm  that  that’s  the  one 
you  want  deleted,  or  ‘8’  to  say  that  its  not  the  one  you  want  deleted.  Leave  the  delete  cycle  by 
keying  ‘8’. 

If  you  wish  to  add  more  to  the  file,  run  it  again,  and  this  time  read  in  the  VRT  file  produced  on 
this  previous  pass. 

At  the  moment,  the  GOD  file  is  the  sum  of  everything  drawn  on  the  screen  during  the  run.  To 
get  just  the  current  description,  you  have  to  pass  the  file  once  more  through  the  program;  making 
no  changes,  just  reading  it  in  and  writing  it  out.  Because  of  this,  I  recommend  NOT  producing  a 
GOD  file  on  the  1“  through  n  —  1"  creation  passes,  saving  it  for  the  last  run,  when  nothing  will 
be  altered. 
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7:  Marimont  Edge  File  Format 

This  section  describes  the  format  of  the  edge  files  created  by  the  Marimont  edge  process. 

7.1  Header 

The  file  is  binary,  and  begins  with  an  128-word  header,  followed  by  a  list  of  linked  edge  lists 
(ledgels);  that  is,  after  the  header,  there  are  a  number  of  items  called  ledgels,  each  of  which  is  a 
list  of  edges. 

The  format  of  the  header  is  as  follows: 

word  interpretation 

0  unused;  ignore 

1  unused;  ignore 

2  bits  per  pixel  in  input  picture 

3  bias  of  laterally  inhibited  picture 

4  “nil-strength”  flag 

5  “end-of-list”  flag 

6  first  row  of  picture  window 

7  first  col  of  picture  window 

8  last  row  of  picture  window 

9  last  col  of  picture  window 

10  1  if  contours  of  positive  regions,  -1  if  those  of  negative  regions. 

11  format  number  of  edge  file  (format  3  currently); 

12-127  unused;  ignore 

7.2  Ledgels 

Next  come  one  or  more  ledgels,  followed  by  a  word  containing  the  end  of  list  flag  (word  5  of  the 
header  above).  Each  ledgel  has  a  two  word  header,  a  list  of  one  or  more  edges,  and  a  word  after 
the  last  edge  containing  the  end  of  list  flag.  Note  that  this  word  flags  the  ends  of  two  kinds  of 
list:  the  list  of  ledgels,  and  the  list  of  edges  that  each  ledgel  contains. 

Each  ledgel  has  the  followi.ig  fo.v.iat: 

word  interpretation 

0  1  if  the  ledgel  is  dosed,  i.e.  the  last  edge  is  linked 

(conceptually)  to  the  first,  else  0. 

1  number  of  left  turns  less  the  number  of  right  turns  in  the  ledgel,  as  one 
traverses  the  ledgel  from  head  to  tail. 

2  beginning  of  this  ledgel’s  list  of  edges 
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[one  or  more  edges] 
last  end'of-list  flag 

7.S  Edgels 

Each  ledgel  has  a  list  of  one  or  more  edges.  There  are  two  kinds  of  edges,  corresponding  to 
“significant"  and  “insignificant”  xero  crossings  in  the  laterally  inhibited  image.  A  significant  xero 
crossing  is  a  transition  of  the  laterally  inhibited  signal  from  positive  to  negative.  An  insignificant 
xero  crossing  is  a  transition  of  the  laterally  inhibited  signal  from  positive  or  negative  (depending 
on  the  sign  of  the  region  whose  encircling  contours  are  being  detected,  see  word  10  of  the  header 
above)  to  xero,  and  as  such  does  not  correspond  to  a  likely  edge  in  the  original  image.  Their 
exclusion  from  linked  edge  lists,  however,  seems  to  destroy  much  useful  connectivity.  So,  we 
include  them  as  “links”  and  assume  that  one  between  two  significant  xero  crossings  is  not  a 
significant  break  in  the  edge.  But  at  insignificant  xero  crossings  we  do  not  estimate  the  strength 
of  the  edge  or  image  intensities  to  cither  side  as  we  do  for  significant  xero  crossings.  Since  less 
information  is  thus  associated  with  insignificant  zero  crossings,  we  use  two  edge  formats  to  save 
space. 

The  format  for  an  edge  is  as  follows: 

word  interpretation 

0  interpolated  row  coordinate  of  edge  point  (floating-point) 

1  interpolated  column  coordinate  of  edge  point  (floating-point) 

2  integer  equal  to  0  if  edge  is  horizontal,  1  if  vertical 

3  two  possible  interpretations:  if  when  interpreted  as  an  integer,  it 
equals  the  nil-strength  flag,  then  this  is  an  insignificant  edge  point, 
the  edge  strength  is  undefined,  and  this  is  the  last  field  in  the  edge. 
Otherwise,  this  edge  is  significant,  this  field  should  be  interpreted  as 
a  floating-point  estimate  of  the  edge  strength,  and  the  next  two  words 
are  fields  belonging  to  this  edge,  all  edge  files  as  of  this  writing  were 
created  by  encircling  positive  regions  (see  word  10  of  the  header),  and 
in  that  case  the  edge  strength  is  always  positive,  and  the  brighter 
side  of  the  edge  is  towards  the  left  as  one  moves  from  the  head  of 
the  ledgel  to  the  tail. 

4  (significant  edges  only!)  estimate  of  original  image  intensity  to  the  left 

5  (significant  edges  only!)  estimate  of  original  image  intensity  to  the  right 
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7.4  Clarifications 

The  edge  location  coordinates  are  in  picture  coordinates,  i.e.  with  respect  to  the  original  picture, 
NOT  w.r.t.  the  window  in  which  they  were  detected.  Integral  locations  correspond  to  pixel 
centers,  not  boundaries,  and  the  topmost,  leftmost  pixel  center  is  taken  to  be  (0, 0),  which  means 
that  an  m  X  n  picture  has  its  upper  left  corner  at  (—0.5, —0.5),  and  its  lower  right  hand  corner 
at  (m  —  0.5,  n  —  0.5)  (since  the  indices  of  the  last  row,  col  are  (m  —  l,n  —  1)). 

For  the  most  part,  edges  are  detected  between  horizontally  or  vertically  adjacent  pixels,  which 
is  why  there  are  only  two  possible  orientations  for  edges.  An  edge  is  detected  between  two  such 
adjacent  pixels  if  the  laterally  inhibited  signal  is  positive  at  one  and  negative  at  the  other.  The 
edge  location  is  estimated  by  the  linearly  interpolated  zero  crossing  between  these  two  values. 
Edge  strength  is  the  difference  between  laterally  inhibited  intensities  one  pixel  to  either  side  of 
estimated  zero  crossing  (in  the  directions  normal  to  that  of  the  zero  crossing);  linear  interpolation 
is  used  to  estimate  each  of  these  intensities,  since  the  points  at  which  they’re  measured  usually 
fall  between  pixels.  The  strength  is  scaled  down  to  correspond  to  the  difference  between  adjacent 
pixels  (i.e.  the  strength  measure  is  laterally  inhibited  contrast  per  pixel  unit  measure).  Image 
intensities  are  estimated  at  these  same  points,  i.e.  one  pixel  to  either  side  of  the  interpolated  zero 
crossing,  in  an  analogous  way,  by  linearly  interpolating  between  each  pair  of  image  intensities. 
However,  both  image  intensities  are  retained. 

The  reason  for  the  “for  the  most  part”  in  the  first  line  of  the  paragraph  above  is  that  “extended” 
zero  crossings  arc  detected  in  the  current  implementation.  An  extended  zero  crossing  is  three 
adjoining  pixels  in  a  row  or  column  with  one  end  pixel  positive,  one  end  pixel  negative,  and  the 
middle  pixel  zero  (values  refer  to  the  laterally  inhibited  signal).  The  estimates  of  edge  location 
and  strength  and  left  and  right  image  intensities  are  made  just  as  in  the  above:  the  location  of  the 
edge  is  an  interpolated  zero-crossing,  the  edge  strength  is  the  (scaled)  difference  in  (interpolated) 
laterally  inhibited  pixel  values  one  pixel  to  either  side  of  the  estimated  location,  and  the  left  and 
right  image  intensities  arc  interpolated  image  intensities  one  pixel  to  either  side  of  the  estimated 
location. 
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8:  Conclusions 

8.1  Summary 

A  principal  research  interest  of  our  group  is  in  developing  a  rule-based  advanced  automated  stereo 
mapping  system  to  function  within  ACRONYM  [Brooks  1981].  Current  mapping  techniques  ignore 
much  of  the  information  available  from  inference  on  single  views  of  a  scene.  This  information 
can  be  useful  for  three-dimensional  surface  interpretation,  and  also  provides  extra  parameters 
for  stereo  matching  (i.e.  surface  orientation,  occlusion  cues).  Our  research  effort  is  directed  at 
establishing  such  monocular  inference  rules  in  a  rule-base  for  stereo  mapping. 

In  deriving  these  rules,  we  perform  analysis  of  both  hand  extracted  and  automatically  produced 
edge  descriptions.  A  facility  has  been  developed  under  the  current  postdoctoral  contract  for 
this  manual  edge  extraction  from  hardcopy  imagery.  We  have  studied  rule  synthesis  for  several 
cases,  including  that  of  orthogonal  trihedral  vertices  -  features  that  dominate  cultural  scenes. 
This  research  is  very  promising,  and  has  shown  the  utility  of  the  rule-based  approach  to  surface 
inference  from  monocular  information.  We  will  be  continuing  our  research  in  rule  synthesis  under 
other  funding. 

Camera  solving  provides  powerful  constraint  on  the  correspondence  problem  in  stereo  matching. 
We  have  developed  a  facility  under  this  contract  for  interactively  registering  images,  determining 
the  parameters  for  transforming  them  (or  their  edge  descriptions)  into  collinear  epipolar  space, 
and  performing  the  actual  image  transformation.  This  determination  is  crucial  to  a  mapping 
process.  Incorporating  an  automated  module  to  provide  data  for  the  camera  solving  is  a  very 
important  next  step. 

We  have  experimented  with  an  existing  stereo  mapping  process,  enhancing  its  flexibility  with 
respect  to  image  format  and  with  respect  to  edge  operator  format,  and  have  been  preparing 
example  outputs  of  its  processing  on  new  imagery.  Our  intent  with  this  effort  has  been  to  show 
the  capabilities  of  a  local  matching  process  and  to  assess  its  applicability  to  the  planned  rule- based 
system. 

Demonatra+ion 

We  will  be  demonstrating  the  following  for  RADC  later  in  May: 

a)  an  interactive  digitizing  test  facility; 

b)  the  use  of  image  edge  descriptions  produced  using  this  facility  and  from 
an  automated  process  [Marimont  1982]  in  synthesizing  rules  for  stereo 
matching; 


c)  crimination  of  Orthogonal  Trihedral  Vertex  (OTV)  inference  rules; 
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d)  development  of  a  system  for  registration  of  image  pairs; 

e)  analysis  of  stereo  imagery  with  an  automated  stereo  process. 

Our  demonstration  will  also  present  related  equipment  and  system  facilities  that 
ol  this  and  other  vision  research. 


we  use  in  support 
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9:  Digitizing  Program 


See  sections  2  and  6  for  documentation  on  using  this  program. 


{$0+}{$H : 20000} 

Program  Tablet( tty , vertf 1 le) ; 


CONST 

CHOP" 100000 ; 

PI  -  3 . 14159265 ; 

DOX  -  432.0:  ODY  •  432.0;  GOO  •  S12.0; 

If*  •  40.0; 

maxdif f erence-0 . 05 :  {  one  twentieth  of  screen  away  } 

{  to  hold  the  lines  associated  with  a  vertex,  incoming  lines  are  negative.  } 
TYPE 

LINESPTR-tlinestype; 

1 inestype-record 

line: integer ; 
nextl ine : 1 Inesptr ; 
end ; 


{  to  hold  the  vertices  with  their  associated  coordinates  and  lines  } 
VERTlISTPTR-fvertliSt; 
vertl  istTecord 

vertex : integer ; 
x ,y : real ; 
header : 1 inesptr ; 
nextvert: vertl istptr; 
end: 

{  to  hold  the  points  associated  with  a  line  } 

POINTPTR-tpoint; 

point*record 

xp ,yp : real ; 
nextpt : pointptr 

end; 

{  to  hold  the  lines  with  their  associated  vertices  and  points  } 
LINELISTPTR-tiinellst; 

1 inel ist*record 

1 1 nenum. in itvert. final  vert: Integer; 
header ipointptr; 
nextl ine: 1  Inal Istptr 
end ; 


TEMPARRAY*array[l . . 20]  of  integer; 

STRNG-packed  array[l:30]  of  ch-.r; 

INTFILE-  FILE  OF  INTEGER;  {  For  graphics  ) 

Var  save*. s avey. XI. Y1.X2,Y2,X3.Y3,X4,Y4, MINX, MINYirs’l ; 

signold . ol Overt , 1 i nenum. pi ch , p icw. vert count , 1 1  recount , verta , vertb : Integer ; 
d*  ,dy ,  gx . gy .  1  s tx .  1  sty .  ttyl  ir.e.boxslze:  integer ;  {  for  drawing  (x,y)  } 

theta. REALX. REALY:real ; 

81,82. answer .dumchar : char ; 

TABLET  10.  Jump  .finished  ,  d(igra,'h .  LAST4  .BOX  FLAG,  boxgod , 

toconfin.r. ::dc;nf Irn.OrC'l'.'X .GOOCP IFX . T,r3ad1n,jdt:'i:tsi?  :sn; 
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vertf lie, inf  lie: text; 

filename: strng; 

vl . startvert , tail wert. 

tptr.vptr:vertl1stptr;  {  vertex  11»t  } 

newlptr  ,  lptr  .  1 , 

scanptr , 1 stscanptr : 1  inesptr  ;  {  lines  of  vertex  list  } 

11 , startline, tailline. 

sptr,savsptr.1stsptr:11ne11stptr;  {  line  list  } 

f  i r stp  , p  , ol dp , ppt r  : pointptr ;  {  points  in  line  list  } 


godfile:  INTFILE;  {  GOD  files  for  plotting  > 

SWITCH:  RECORD  {  Little  hack  for  putting  reals  Into  integers  } 

CASE  F 00 : 800LE AN  OF  {  for  60D  files  } 

TRUE:  ( I : INTEGER) ; 

FALSE:  (R: REAL) 

END; 

{  external  DD  graphics  procedures  in  GPAS . REL[TST , AAM]  } 
procedure  ginit:  extern;  {  init  } 

procedure  scrset;  extern;  {  clear  screen  } 

procedure  width( w: Integer ) ;  extern;  {  no-op  } 
procedure  drken;  extern;  {  make  line  dark  } 

procedure  liten;  extern;  {  make  line  light  } 

procedure  inven;  extern;  {  swap  above  } 

procedure  move( x ,y : integer ) ;  extern:  {  invisible  line  to  x,y  } 

procedure  draw( x , y : integer ) ;  extern;  {  visible  line  to  x.y  } 

procedure  dpyup;  extern;  {  draw  to  screen  } 

procedure  inchar(var  key , ctl key .metakey : integer) ;  extern; 
procedure  b  i  tin i ( 1  ine :  integer ) ;  extern;  {  init  tty  for  input  } 

procedure  bitpad(var  resul 1 . resul 2 . resul3 : integer ) ;  extern;  {  tablet  input  } 

function  ttytab(var  ;sul 1 , resul 2 , resu!3 : Integer ): boolean ;  extern;  {tablet  or  tty  test  ) 

procedure  ppset;  extern;  {  set  page  printer  } 

procedure  ppdone;  extern;  {  release  page  printer  } 


{  REAOTABLET  DISTANCE  INITTEMPARRAY  GRAPHICS  (DD  AND  GOD)  } 

Procedure  REAOTABLET( var  B:char;var  X,Y:real); 
var  T , IX . IY : integer ; TA : real  ; 

{  This  procedure  reads  jn  data  from  either  the  terminal  or  the  tablet. 

The  data  consists  of  which  button  was  pressed,  and  the  coordinates 
of  the  point  at  which  it  was  pressed.  The  boolean  TABLETIO  Indicates 
whether  the  procedure  is  reading  from  the  tablet  or  from  the  terminal.  } 

Begin 

if  TABLETIO  than 
begin 

bitpad(T.IX.IY); 

{  IF  TTYTAB(T.IX.IY)  THEN  WRITELN( TTY. 'tty  activity’);} 

{  wr 1  tel n( tty . t :  0 .  *  '.1x:0,'  *,1y:0);  } 

B:-chr(4S+T) : X : -IX; Y: -IY; 
end 

else 

begin 

readln(tty) : 

read( tty ,B) ; read( tty , X) ;read( tty ,  Y) ; 
end; 

End; 

PROCEDURE  GETCOORDS; 

Begin 

READTABLET(31,*l.Yi); 
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REALX  :•  ( Xl-MINX )/p lew;  REALY  :•  ( Yl -MINY )/piCh ; 

IF  REALX  <  0.0  THEN  REALX:-0.0 

ELSE  IF  REALX  >  1.0  THEN  RE'ALX:*1.0; 

IF  REALY  <  0.0  THEN  R£ALY:*0.0 

ELSE  IF  REALY  >  1.0  THEN  R£ALY:-1.0; 

{  write ln(tty. -B1*  * ,B1,  -,  X1*’,X1:0,'(* .REALX : 5 ; 3 , ’ ) .  Yl--.  Y1;0 

End; 


FUNCTION  ASK(N: integer) :  BOOLEAN; 

VAR  goodans : BOOLEAN ; 

Begin  {  ASK  } 
goodans : ‘FALSE ; 

Repeat 

case  N  of 

1 : WRI TE ( TTY , 'Want  Tablet  10  (y/n)?:-); 

2 : WRITE ( TTY , 'Want  graphics  output  (y/n)?:*); 

3 : WR I TE ( TTY . '  to  the  00  {y/n)?;*); 

4 : WR ITE ( TTY ,  '  to  the  GOOFILE  (y/n)?:1); 

5 : WRITE ( TTY , 'Want  to  read  in  a  file  (y/n)?:-); 

6 :WRITE(TTY, 'Happy  with  this  frame  (y/n)?:-); 

7 : WRITE( TTY . ‘ Want  to  DELETE  segments  (y/n)?:-); 

8 : WR I TE ( TTY , 'Do  another  deletion  (y/n)?:-); 

9 : WRITE( TTY .  '  I s  this  the  line  to  be  deleted  (y/n)?:-); 

10  ;WRITE( TTY , '  Want  boxes  draw  around  vertices  (y/n)?:-); 

11 : WRI TE( TTY , '  want  frame  drawn  in  GOD  file  (y/n)?:-); 

12 : WRI TE ( TTY , ’  Want  to  add  to  file  from  TABLET  ( y/n ) ? : - ) ; 

End  ; 

readln( tty) ; read( tty . answer ) ;  if(answer  -  -y-)  or  (answer  •  -Y- 
Begin  ASK TRUE ; goodans : -TRUE ;  End 
else  Begin  ASK : -FALSE ; goodans : -TRUE ;  End; 

UNTIL  goodans; 

End:  {  ASK  } 

Function  DISTANCE(X1,  Yl.  X2 ,  Y2:  Real):Real;  {Finds  the  distance  bet. 
begin 

distance:-  Sqrt  (Sqr  (X2  -  XI)  +  Sqr  (Y2  -  Yl)); 


Procedure  INI TTEMPARRAY( var  a : temparray )  ; 
var  i : integer  ; 
begin 

for  i:-  1  to  20  do 
a[ 1 ) : -0 ; 

end ; 

Function  HAX( a , b : Integer ) :  INTEGER; 

Begin 

IF  A  >  B  THEN  MAX  :■  A  ELSE  MAX  :•  B; 

End; 

{  Routines  to  make  GOD  files  } 

PROCEDURE  0UT(N: INTEGER) ; 

BEGIN 

GODFILEt  :-  N:  PUT(GOOFILE) 

END; 


PROCEDURE  OUTR£AL(V:R£AL); 

BEGIN 

SWITCH. R  :•  V; 

OUT (SWITCH . I ) 

END; 

PROCEDURE  PLOTLINE (XI. Y1.X2.Y2: REAL ) ; 
BEGIN 

OUT (  7 )  ;  {  LI»'E  } 

OUTREAL(Xl):  OUTCr'-f  Yl ) ; 


, -(-.REALY:5:3,-)-);  ) 


)  then 


2  pts.  } 
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OUTRE AL( X2 ) ;  OUTREAL(Y2); 

OUT { 0 )  {  THICK  > 

END; 

{  Procedures  to  display  graphics  on  the  OD  terminal.  } 


Procedure  BOX; 

Begin  {  BOX  } 

IF  BOXFLAG  THEN 
Begin 

IF  ODGRAFX  THEN 
Begin 

MOVE( dx-boxslze,dy-boxslze); draw) dx-boxsize, dy+ box  size) 
0RAW(dx+boxsize.dy+boxsize);Draw(dx+boxs1ze,dy-boxs1ze) 
Dr a w( dx-boxsize, dy-box size) .MOVE) dx . dy ) ; 

IF  NOT { RE ADIN )  THEN  OPYUP ; 

End; 

IF  GOOGRAFX  THEN 
Begin 

P LOTLINE) GX-boxsize. GY -boxsize.GX-boxsize.GY+boxsIze); 
PLOTLINE(6X-boxsize.GY*boxsize,GX+boxs1ze.GY+boxs1ze): 
PLOTLINE(GX+boxsize.GY+boxsize,GX+boxsize,GY-boxs1ze); 
PLOTLINE(GX+boxsize,GY-boxsize,GX-boxsize.GY-boxsize); 

End ; 

End ; 

End;  {  /OX  } 

procedure  DARKDRAW; 

Begin 

IF  GODGRAFX  THEN 
Begin 

gx:=round( realx*GOD)  : 
gy : '  round) realy*GOO) : 

End ; 

IF  DDGRAFX  THEN 
Begin 

dx  :  -  round) real x*DDX  +  l f x) ; 
dy  ;•  round) ( 1 . O-realy )*DDY) ; 

MOVE ( dx . dy ) ; 

IF  NOT ( REAOIN )  THEN  OPYUP; 

End;  {  move  if  DO,  nothing  if  GOD  } 

End ; 

procedure  LIGHTDRAW; 

Begin 

IF  GOOGRAFX  THEN 
Begin 

lstx:*gx;  1  sty : - gy ; 
gx:  *round(real x’GOD) ; 
gy : *  round) realy*GOD) ; 

PLOTLINE (Istx.lsty.gx.gy); 

End; 

IF  DOGRAFX  THEN 
Bagin 

dx  ;  *  round) real x*OOX+lf x) ; 
dy  :•  round) ( 1 . O-realy )*ODY) ; 

ORAW(dx.dy); 

IF  NOT (REAOIN )  THEN  OPYUP; 

End; 

End ; 

PROCEDURE  FR.V'LR/J: 

Begin  {  FRAlxc  } 

realx:-0.0;realy:-0.0;OA.TKOR,\W;  {  DRAW  FRAME  AROUND  SCREEN  ) 
realx:-0.0:realy:*1.0 ; LIGHTDR  ’  W; 
rea1x:*1.0;rea1y:*1.0 ;  LIGHfOR.-  1: 
real  x :  -1 . 0 ;  re  Jly :  *0 . 0 ;  LIC’ITLh."  I; 

-•llr.:-O.C;rdJly:,-O.C.M'‘  •*:’  ’; 
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End ; {  FRAME  } 

PROCEDURE  GOOINI; 

Begin  {  GODINI  } 

R£WRITE(GODFILE , 'TAB. GOO'); 

OUT ( 97 )  ;  {  DDINIT  } 

OUT ( 1 ) ;  {  SCREEN  ) 

OUTRE AL (O.O); 

OUTRE AL(O.O); 

OUTREAL(GOO) ; 

OUTREAL(GDO) ; 

OUT ( 4 ) ;  {  LITEN  } 

BOXGOD : »  ASK( 11); 

If  BOXGOD  THEN  f RAMEDRAM ; 

End.  {  GODINI  ) 

PROCEDURE  GODFINI: 

VAR  I : integer ; 

Begin  {  DPYUP  and  KILJ08  } 

OUT ( 17)  ;  {  OPYUP  } 

OUT(-l)  ; 

OUT ( - 1 ) ; 

fOR  I  1  TO  32  DO  OUT(O)  {  KILJOB  } 
End;  {  DPYUP  and  KILJOB  } 


PROCEDURE  DDINI; 

Begin  {  DDINI  ) 
PPSET; 
ginlt : 
scrset ; 

ORKEN; 

IF  not ( GODGRAF X ) 
End;  {  DDINI  } 


{  set  up  small  page  printer  for  text  } 


{  dark  here  so's  works  like  SUN  } 
THEN  FRAMEORAW;  {  frame  if  Just  DD  } 


PROCEDURE  DOF  INI ; 

Begin  {  page  printer  clear  } 
PPDCNE  ; 

End;  {  page  printer  clear  ) 


{  INITLINEIIST.  INITVERTLIST .  ENTERLINE,  ENTERVERTEX,  ADDLINEPTR  } 
Procedure  INI TLINELIST  ; 

(this  initializes  the  list  of  lines) 

Beg  in 

new(ll); 

1 1  * . nextl ine: -nil; 
lit. header  : »ni  1  ; 
s tar  1 1 ine : • 1 1  ; 
ta i 1 1 i ne :  *1 1 ; 
fnd;{INITLINELIST) 


Procedure  INITVERTLIST; 

(this  initializes  the  list  of  vertices) 
Begin 

new( v  1 ) ; 

v 1 1 . nextver t : «n  1 1 ; 
v l * . header  :  - n  1 1 ; 
startvert : *vl ; 
tai 1  vert : »vl  ; 

End;{INITVERTLI$T) 


Procedure  EN1 ERLI  E  ; 

(this  cuts  the  line  Into  li.-.olidt) 
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Begin 

linecount:-linecount+l: 
tailline*. linenum.-linecount; 
tai 1 1 iner . initvert:-verte; 
tail  1  met .  f  malvert:-vertb; 
tai  1 1  inet  header:-f irstp; 
f  irstp : *ni 1 ; 
new(  11); 

taillinet.nextline:*ll; 
tai  1 1  me  :  -1 1  ; 
taillmet.ne*tline:*nil; 
tnd;  {ENTERLINE} 


Procedure  ENTERVERTEX; 

{  starts  a  new  vertex  } 

Beg  l  n 

tai 1 vertt . % : -REALX ; 
ta l 1 ver t t . y : -REALY ; 
vertcount:*vertcount+l; 
tailvertt.vertex:-vertcount; 
ta 1 1 ver 1 1  header  . -nil  ; 
new( v 1  ) ; 

tailvertT.nextvert:*vl  : 
ta i 1 ver  t : *vl ; 
t a i 1 ver  t ' ,nextvert:*n11  ; 

End;  {  ENTERVERTEX  } 


Procedure  ADDLINEPTR(1 i ne . ver tnum : In teger )  ; 

{  this  adds  a  line  to  the  linked  list  of  lines  associated  with  a  vertex  } 
Begin 

vl ; - s  tartvert : 

while  vl*. vertex  <>  vertnum  do  vl : -vl t . nextvert ; 

1 ptr : *  vl ♦ . header  ; 
if  1 p  tr  =  nil  then 
begin 

new( newl ptr ) ; 
vl*.header:»newlptr; 

end 

else 

begin 

while  1 ptrt . nextl  ine  <>  nil  do  lptr : -Iptrt . nextl Ine; 
new( newl ptr ) ; 

Iptrt. nextline: -newl ptr ; 

end ; 

newlptrt . l ine ; -1 ine: 
new Iptrt. nextl 1ne:-nil  ; 

End;  {ADOLINEPTR} 


Procedure  F INDOLDVERTEXf var  vp tr : ver 1 1  is tp tr ; var  ol dvert : integer ; xol d ,yol d: real ) ; 

{  This  fines  the  old  vertex  which  is  nearest  to  the  point  at  which  BZ  was 
pressed  If  the  nearest  old  vertex  is  further  than  maxdlfference  away,  It 
tells  how  far  away  it  is.  } 

var  vertnum: integer ; 
tptr : ver tl 1 stptr ; 

Begin 

tptr : -startvert; 

vptr : -startvert ; 

while  ( tptrotai  1  vert)  do 

begin 

1 f  DISTANCE ( tptrt . x, tptrt .y , xol d , yol d)<DIST VICE (vptr t .x, vptrt .j , xold.ypl  1 


vptr:*tptr; 
tptr : -tptrt . nextvert; 

end : 

old wart : -vptrt . vertex ; 

if  (OISTANCE(vptr r.x, vptrt. y . xol d ,yold)>maxdiff erence)  then 

writeln(tty, 'the  nearest  old  vertex  is  * ,DISTANCE( vptrt . x . vptrt .y.xol d ,yold 

2  ,  '  away ' ) ; 

End;  {FINOOLDVERTEX} 

PROCEDURE  AIMVERT(vert: integer) ; 

Begin 

tptr:*startvert; 

while  tptrr  .  vertexovert  do 

tptr:-tptrf. nextvert; 

R£ALX:-TPTRr.X; REALY : -TPTRf . Y ; 

End ; 

FUNCTION  F IND0LDP0INT( var  re tsptr ; 1 i nel 1 S tptr ; TESTX . TESTY : real ) :  BOOLEAN; 

{  This  finds  the  old  point  which  is  nearest  to  the  point  at  which  B2  was 
pressed.  If  the  nearest  old  point  is  further  than  maxdif ference  away,  it 
tells  how  far  away  it  is.  } 

VAR  OLDOIST : REAL ; 


PROCEDURE  SEARCHPOINT(  vl  ptr  ;lir,esptr;xold.yold:  real  )  ; 

Var  tstsegnum; integer ; 

Begin 

tstsegnum: -ABS(vlptrt .  1  ine) ;  {  get  linenumber  } 

sptr : -startl ine; 

while  sptrt .  1  inenumotstseqnum  do 
sptr:*sptrt.nextl1ne; 
pptr  :  -sptr* . header  ; 

WHILE  pptroNIL  00 
Begin 

IF  DISTANCE( pptr t , xp , pp trt ,yp , xold .yol d )  <  OLDDIST  THEN 
Begin 

OLDDIST  :*  DISTANCE ( pp trr . xp , pptr t . yp , xol d . yol d ) ; 
savsptr : -sptr ; 

End; 

pptr : -pptrr . nextpt ; 

End; 

End ; 


Begin 

AIMVERT(verta) ; 
savsptr ; *ni  1  ; 
scanptr : -tptrt . header ; 

OLDDIST:* 1.0; 

REPEAT 

SEARCHP0INT( scanptr .TESTX, TESTY) ; 
scanptr  :  -  scarptr-r .  next!  Ine; 

UNTIL  scanptr  •  nil  ; 

IF  OLDOIST  >  maxdifference  THEN 
Begin 

{  WR I TELN( TTY ,  ’No  such  point  on  a  line  from  that  vertex');  } 

FIHD0LDP01NT: -FALSE; 

End 

ELSE 

Digin 

verta:-  savsptr*.  ini  tv  trt: 
vertb:-savsptr*.f1n3lv:rt; 

retsptr : -savsptr ;  {  copy  bsck  for  return  } 

F  INOOLCrO.  r  :  - T R ‘J E  ; 

End: 

End;  {FI  .CCL  . -0I"T) 


PPCCEO'.'OE  r.  -  EST(s3^ptr:Ht?li  ;tdtr;pr;cffdur3  C".'.?MIC;  doi  t :  bool  aen) ; 
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Begin 

verta: -segptrt . Ini tvert ; 
vertb ■ -segptrt . final  vert ; 

GRAPHIC; 

AIMVERT(verta)  ;  DARKDRAW;  {  start  drawing  } 
pptr : -segptrt . header ; RE AO  IN  :  -  TRUE ; 
while  pptr  <>  nil  00 
Begin 

REALX  :  -  pptrt . xp ; REALY : -pptrt  .  yp; 
pptr : -pptrt . nextpt ; 

LIGHTORAW; 

End; 

AIMVERT(vertb) ; LIGHTORAW ; REAOIN : -FALSE  ; 

DPVUP ; DRKEN ;  {  finishing  drawing  } 

END; 

FUNCTION  CHECKOLDVERTEX(var  r e t sp t r : 1  i  ne 1  i  s t ptr  ) :  BOOLEAN; 

Var  tstsegnum. checksegnum: integer; 
tstvp  tr : vert  1  is tptr  ; 
scnvl ptr ; 1 inesptr  ; 

TSTT  BOOLEAN: 

FUNCTION  CHECXVERTSOFLINE:  BOOLEAN; 

Begin 

tstsegnum : *A8S( scanptrt  .  1  ine)  :  {  get  linenumber  } 

scnvlptr:*tstvptr’. header;  {  scan  vertex  line  list  } 

checksegnum : -  -  1 ; 

WHILE  ((scnvlptr  <>  NIL )  ANO  (tstsegnum  check segnum) ) 
DO  Begin 

checksegnum.-ABSfscnvlptrt.line); 

scnvlptr;-scnvlptrt.nextline; 

End ; 

IF  tstsegnum  <>  checksegnum  THEN 
Begin 

CHECXVERTSOFLINE : -FALSE ; 
scanptr:=scanptrt. next  line; 

End 

ELSE  CHECXVERTSOFLINE : -TRUE ; 

End ; 

Begin 

retsptr  :  »n  i  1  ; 

F INDOLOVER TE  X ( t s tvp tr , ver tb , RE ALX , RE ALY ) ; 
sav sptr : =n  i  1 ;  {  tptr  and  tstvptr  } 

scanptr:-tptrt. header; 

TSTT: -FALSE; 

WHILE  ( ( scanptr  <>  NIL)  ANO  NOT(TSTT)) 

00  TSTT; -CHECXVERTSOFLINE; 

IF  TSTT  THEN 
8eg1n 

sptr : -start! ine; 

while  sptrt .  l  inenumotstsegnum  do 
sptr  :  -sptrt . nextl Ine; 
retuptr : -sptr ; 

CHECXOLDVERTEX  :•  TRUE; 

End 

ELSE  CHECXOLDVERTEX  :■  FALSE 
End;  {  CHECXOLDVERTEX  } 

PROCEDURE  DELETESEG( segptr : 1  Inal Istptr) ; 

PROCEOURE  OELSFROMV( vertnum: Integer); 

Begin 

A INVERT ( ver  tnum) ; 
scanptr : -tptrt . header ; 
lstscanptr: -tptr t , header ; 
while  ABS( scanptrt .  1  i,ie)  <>  11, mum  CO 
Bagin 

1  stscciptr  ;  :tr  ; 


scanptr:-scanptrt. nextl In#; 

End; 

IF  scanptr-tptrt . header  than  tptrt.headar:-seanptrt.nextl1ne 
ELSE  Istscanptrt .nextl Ine: -scanptrt .naxtllne; 

End; 

Begin 

verta;  -segptrt .  inltvert; 
vertb ; •segptrt. final  vert; 

1 inenum: -segptrt . l inenum; 

OELSFROMV(verta) ;  {  removal  Ine  from  start  vertex  Hat} 

OELSFROMV(vertb) ;  {  removal Ine  from  end  vertex  list) 

sptr:-startl1ne; 

1 stsptr : -nil ; 
while  sptr  <>  segptr  00 
Begin 

1 stsptr : -sptr ; 
sptr:*sptrt. nextl Ine; 

End ; 

IF  lstsptr-nll  THEM  startl ine : -sptrt . nextl Ine 
ELSE  1stsptrt.nextl1ne;-sptrt. naxtllne; 

End;  {  DELETESEG  } 

PROCEDURE  DELETEPHASE; 
var  segptr : 1 Inel Istptr ; 

PROCEDURE  SETDELETION; 

Begin 

F INDOLDVERTEX ( tptr, verta, REALX , REALY) ; 

GETCOORDS;  {  select  line  or  vertex  point  } 

T : “FALSE ; 

CASE  BT  of 

•1':  T:-  CHECKOLDVERTEX( segptr ) ; 

’  4 ' :  T  :  •  F INDOLDPOINT ( segptr , RE ALX , REALY) ; 
others:  Wr 1  tel n( tty .  '  1  for  vertex;  4  for  point') 

End ; 

IF  T  THEN 
Begin 

REM0VETES1 (segptr. LI  TEN. TRUE); 

NEEOCONF IRM : *TRUE ; 

wr1teln( tty , ’<B*2>  for  confirm,  <8-8>  for  reject  deletion’); 

End 

ELSE  wr 1  tel n( TTY , 'Coul dn ' ' t  delete  line  (from  vertex  verta : 0 .')’) ; 

End;  {  SETDELETION  } 

PROCEDURE  ERRMSG; 

Begin 

wr1teln(TTY, 'B»l  or  4  for  deletion.  B-Z  for  confirm,  B*B  for  reject/quit’); 

End; 

Begin  {  Deletes  segments  } 

PPSET;  ' 

wr 1  tel n( TTY . ’ Sel ect  vertex  <6-l>,  then  vertex  <B*1>  or  point  on  line  <B-4>;  <B-8>  quits.') 
wr  1  te( TTY ,  '  Once  selected,  <B-Z>  confirms  deletion.  <B*B>  cancels  It.’); 

FINISHED: -FALSE; NEEOCONF IRM; -FALSE; 

REPEAT 

GETCOORDS;  {  select  vertex  } 

CASE  B1  OF 

'1':  IF  NOT (NEEOCONF IRM)  THEN  SETDELETION 

ELSE  ERRMSG;  {  CHOSE  THE  LINE  TO  DELETE  } 

'2':  IF  NEEDCONFIRM  THEN  {  THIS  IS  CONFIRMATION  } 

Begin 

OELETESEG( segptr ); NEEOCONF IRM: -FALSE; 

wr1teln(tty,  'For  deletions  B-l  then  1  or  4:  b-8  Is  quit'); 

End 

ELSE  ERRMSG; 

'8':  IF  NEEDCONFIRM  THEN  {  THIS  IS  REJECTING  THE  DELETION  } 

Begin 

REMOVE  TEST(  segptr.  DRKEN, FALSE )  ;KSEDC0"FIRf1:  -r.'.'.SS; 
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wr1te1n( tty , *  For  deletions  B- 1  then  1  or  4;  B*8  is  quit') 
End 

ELSE  F inished : -true ; 

OTHERS:  ERRMSG 

END;  {  OF  CASE  ) 

UNTIL  FINISHED; 

T  : -TRUE  FINISHED: -FALSE  ;LAST4: -FALSE; 

»r i tel n( tty Return ing  to  entry  mode’); 

End:{  Deletes  segments  ) 

PROCEDURE  LASTDELETE;  {  last  chance  at  deleting  } 

Begin 
PPSET ; 

Writeln(TTY. 'DELETE?  <8*2>  for  YES.  <8-8>  for  NO'); 

T: -FALSE; 

REPEAT 

READTABLET (B1 .  X4.  Y4); 

CASE  B1  of 

•2':  OELETEPHASE ;  {  do  deletions  ) 

•8’ :  T : -TRUE ; 

others:  wr i tel n( TTY . ' <B-2>  for  deletions  or  <B-6>  for  none*) 

End ; 

UNTIL  T; 

End ; 


Procedure  STARTOLOVERTEX ; 

{  This  procedure  finds  an  old  vertex,  enters  the  number  of  the  incoming 
line  in  the  linked  list  of  lines  associated  with  it.  enters  the  number 
of  the  old  vertex  in  tne  line  buffer,  and  calls  ENTERLINE  to  move  the 
contents  of  the  buffer  into  the  linked  list  of  lines.  It  then  draws  e 
line  to  the  old  vertex.  ) 

var  1  inenum , ol dvert : in teger ; 

Begin 

toconfirm:*false:LAST4:-FALSE; 
if  jump  then 
begin  {  jump  } 

F I NDOLDVERTEX ( vptr.oldvert, RE ALX , REALY) ; 

REALX:-vptrt.x:REALY:-vptr».y; 

verta : »ol dvert ; 

jump : *f  al se : 

if  ddgraph  then  DARKDRAW ; 
end  {  Jump  } 

else 

begin  {  no  jump  } 

F INDOLDVERTEX( vptr, ol dvert, REALX, REALY); 
ver tb : -ol dvnrt ; 

REALX:-vptrt.x ; REALY : -vptrt ,y ; 

ENTERLINE; 

ADOLINEPTR( linecount, verta); 

ADDLINEPTR( -1 inecoun t , ver tb ) : 
if  ddgraph  then  LIGHTDRAW; 
verta : -ol dvert ; 
end;  {  no  Jump  } 

End;  {STARTOLD. TRTEX} 


Procedure  STAR.  .  £.:7t.X; 

{  This  p.-jr:.'  :  cr:  tpp  i  :.;:i  pitry  in  linked  list  of  v;.‘t1c.s 
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which  contains  tha  number  and  coordinates  of  the  vertex,  and  -Jds 
the  vertex  number  to  the  line  buffer.  If  the  vertex  Is  at  the  end 
of  the  line  (if  84  has  bean  previously  pressed  to  add  points  to  that 
line).  It  calls  ENTERLINE  to  add  the  contents  of  the  line  buffer 
to  the  linked  list  of  lines.  It  also  draws  a  line  to  the  vertex.  } 

var  1 inenum: Integer: 

Begin;  {STARTNEWVERTEX} 

toconf 1 rm;  *f  al se ; LAST* ; -FALSE ; 
if  jump  then 
begin  {  jump  } 

ENTERVERTEX; 
verta: -vertcount: 
jump: -false; 

if  ddgraph  then  Begin  DARKDRAW.-BOX  End; 
end  {  jump  } 
else 

begin  {  no  Jump  } 

ENTERVERTEX; 
vertb : -vertcount; 

ENTERLINE; 

ADDLINEPTR( 1 inecount , verta ) ; 

AOOLINEPTR( - 1 inecount. vertb); 
if  ddgraph  then  Begin  LIGHTORAW-.BOX  End; 
verta: -vertcount; 
end;  {  no  jump  } 

End;  {STARTNEWVERTEX} 


Procedure  NEWPOINT; 

{  This  procedure  adds  the  coordinates  of  a  new  point  to  the  linked 
list  of  points  which  are  associated  with  a  vertex.  } 

Begin 

IF  toconfirm  then  OELETEPHASE 
ELSE 

IF  NOT ( JUMP )  THEN 
Begin 

If  firstp-nll  then 
begin 
new(p); 
f irstp:-p; 
oldp:-p 
end 

else  begin 

new(p) ; 

oldpt.nextpt:-p; 
ol dp : *p ; 
end ; 

pt.xp:*REALX;pf.yp: -RSALY ; 
pt .nextpt:*n11 ; 

If  ddgraph  then  LIGriTDRA'J; 

LAST*: -TRUE; 

End 

ELSE  WRITELN{ TTY , ’ Put  the  pen  dorm  first’); 
toconf 1rm:*falte; 

End;  {NEWPOINT} 


{  UTILITIES  } 

FUNCTION  REAOIHF ILE :  “»0‘  '."'.X; 

VAR  basel  1ne,bas:v:rt ,  I , m..'.v:rtj,  •.  •  T,n-1  Ir.:; .  1  vr.i.ol  J1 1  r.e , QV.v:T :  I.’TEGER; 
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TVAL: BOOLEAN; 
tptr ; vertl Istptr ; 

Procedure  READVERTS; 

{  This  procedure  reads  the  VERTICES  from  the  vertflle  } 

Begin 

basel Ine: *  1 Inecount ; basevert : -vert count ; 
read(infi1e. dumchar , dumchar , dumchar ) ; 
read( inf ile. numverts) ;  {  get  past  'V  •'  for  N  } 
read! n( Inf  1 1  a) : 

FOR  I;*  1  to  numverts  00 
Begin 

read( inf  11 e , vnum, dumchar .dumchar . REALX . dumchar .REALY , dumchar , dumchar) ; 

{  number,  coorda,  coordb  } 

ENfERVERTEX; 

while  not(eoln(1nf11e))  DO 
Begin 

re ad {Infile. oldllne); 

IF  oldline  <  0  then  s1gnold:»-l  else  s1gnold:*l; 

oldl ine ; *ABS( oldl Ine) ;  oldl Ine: •( ol dl Ine+basel 1ne)*s1gnold; 

ADDLIN£PTR(oldl ine , vnum+bas evert ) ; 

End ; 

read1n(  inf  He) ; 

End ; 

vertcount  :•  basevert  +  numverts; 
linecount  :«  baseline; 

End;  {  REAOVERTS  } 

Procedure  READLINES; 

{  This  procedure  reads  the  LINES  from  the  Inflle  } 

PROCEDURE  FINDANDOR AW( ver t : 1 nteger ; procedure  FUNC); 

Begin 

tptr.-  -startvert: 
while  tptrt .  vertexovert  do 
tptr  :* tptr t . nextvert ; 

REALX: >TPTRr . X ; RE ALY : -TPTRt . Y ; 

FUNC: BOX; 

End ; 


Begin 

readln( inf ile, dumchar .dumchar. dumchar. numl Ines) ;  {  get  past  'V  •*  for  N  ) 
FOR  I:-  1  to  numl ines  00 
Begin 

read(infile. 1 n urn, dumchar , dumchar . verta , dumchar , ver tb .dumchar, dumchar ) ; 
ver la ; -ver tarda  sever t ; vertb . -ver tbrbas evert ; 

IF  ddgraph  THEN  F INDANODRAW( verta .DARXDRAW) ; 
f  irstp  :  -ni  1  ; 

while  not(eoln(1nf lie))  DO 

Begin 

read( inf i le, dumchar , REALX . dumchar .REALY, dumchar) ; 

NEWPOINT; 

End ; 

IF  ddgraph  THEN  F IN0AND0RAW( ver tb , LIGHT DRAW) ; 

ENTERLINE; 
readln( inf ile) ; 

End ; 

End; {HEADLINES} 


Begin 

TVAL ; -ASK( S) ; 

IF  TVAL  THEN 
Begin 
PPSET ; 

{  WRITE(TTY. •Filename;’);  } 

{  READLN( TTY . f ilenane : QUANT ; [  ]);  } 

{  WRITELN{ TTY .QUANT : 0) ;  } 

RESET( Inf  lie , ’NEU. VYT ' ) ; 


readin:-true; 

REAOVERTS: 

HEADLINES; 

IE  (ddgraph  and  ddgrafx)  THEN  OPYUP; 
End ; 

R£ADINFILE:-TVAL; 

End;  {  REAOINFILE  } 


Function  MAG( XI .  Yl.  XZ,  Y2:  Real ): INTEGER ;  {  Finds  the  distance  bet.  2  pts.  } 

BEGIN 

MAG:-  round(Sqrt  (Sqr  (X2  -  XI)  ♦  Sqr  (Y2  -  Yl))); 

END; 

Function  ANGLE  (XI.  Yl.  X2.  Y2;  Real):  Real;  {  Finds  angle  theta  bet.  2  lines  } 
VAR  x.Y, xlt -.real ; 

BEGIN 

X  :-  X2  -  XI; 

Y  :•  Y2  -  Yl; 

IF  abs(X)  >  abs(y )  THEN 
xLT  :»  Arctan  (Y/X) 

ELSE  xLT  :-  (PI/2)  -  Arctan  (X/Y); 
angle  : ■  xl t ; 

END: 

FUNCTION  INITTABLET : BOOLEAN; 

Begin 

Wrlteln  ( tty ,' Lower- 1 ef t  corner;  '): 

REAOTABLET (B1 ,  XI,  Yl); 

Wrlteln  (tty.  'Top-left  corner;  '); 

REAOTABLET ( B 1 ,  X2 ,  Y2  ) ; 

Wrlteln  (tty,  'Top-right  corner;  '); 

REAOTABLET ( B 1  .  X3.  Y3): 

Wrlteln  (tty,  'Lower-right  corner;  * ) ; 

REAOTABLET ( B1 ,  X4.  Y4); 

{  'nese  calculations  are  dumb  ) 
minx  : *  XI :  miny  : *  Yl ; 

Theta  :-  ( ( ANGLE ( X 1 .  Yl ,  X2,  Y2)  -  Pi/2.0)  +  ANGLE(X2 ,  Y2.  X3,  Y3)  ♦ 

( ANGLE ( X4 .  Y 4 ,  X3.  Y3)  -  Pi/2.0)  +  ANGLE(X1 ,  Yl,  X4,  Y4))  /4.0; 

pich  ;*  MAX ( MAG ( XI ,  Yl,  X2,  Y2),MAG(X4,  Y4,  X3,  Y3)); 
picw  :*  MAX ( MAG( X 1 ,  Yl,  X4.  Y4).MAG(X2,  Y2 ,  X3,  Y3 ) ) ; 

Wr  i  te( tty  .  'Theta-' . theta : 5 : 3 , '  height-'  ,pich:0,  '  width- ’ , picw: 0) ; 
Writeln(tty.  '  <B*2»  confirms,  <B-8>  rejects'); 

REAOTABLET ( B1 ,  X4.  Y4); 

CASE  B1  of 

'2':  INITTABLET: -TRUE ;  {  make  sure  happy  with  this  } 

'8'  INITTABLET:-FALSE; 

others:  wr 1  tel n( TTY . ' 2 ( conf Irm)  or  B(reject)') 

End ; 

End ; 

PROCEDURE  WRITEFILE ; 

TYPE  bigarray-arrey[l. .500]  of  integer; 

VAR  VTRAN.LTRAN.blgarrey; 

VCOUNT , LCOUNT : integer; 

PROCEDURE  BUILDVLTASLE ; 

Begin 

tptr:-stirt/ort; 

VCOUNT: -0: 

while  ( totrotell vert)  do 
Bsgin 

verta : -tptr* . vertex ; 

IF  tptrt. header  <>  nil  then 
Begin 

’/COUNT :  -VCOUM+l ; 

VTRAH[verta]:-VCOJiT; 

End 

VT 1 '  t'[  v  tr  1 1] :  -  - 1 ; 


{  > 
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tptr:*tptrf.nextvert; 

End; 

sptr:>startl1ne; 

LCOIINT  :*0; 

while  sptr  <>  tallllne  do  {  no  nulls  In  saga.,  they've  baan  removed  } 

Begin 

1 Inenum: *sptrt . 1 Inenum; 

LCOUNT ; -LCOUNT+l ; 

LTRAN[1 Inanum] : -LCOUNT 
sptr : -spirt .next 1 ine; 

End; 

End; 


Begin  {  WRITE  OUT  VERTICES  AND  EOGES  } 

{  write  file  with  entries; 

V  *  the  number  of  vertices, 
cvertex  entries>  (N  of  them), 

L  ■  number  of  lines, 

<line  entr1es>  (M  of  them).  } 

BUILOVLTABLE ;  {  make  index  transforms  for  vertices  and  lines  } 

wr i tel n( tty . ’ V  *  1 , VCOUNT ; 0 . • .  L  *  '.LCOUNT:0); 

{  do  Vertices  } 

wr i teln( vertf i 1 e , ’ V  *  ',VC0UNT:0,'  height  •  '.p1ch:0,',  width  •  ',p1ew:0,’  (mm)1); 

vl : -startvert ; 

repeat 

verta ; =VTRAN[vl t .vertex] ; 

IF  verta  >0  then  wrl te( vertf lie, verta: 0 ,vl* . x :0 : 3 ,vlt .y : 0 : 3 ; 

1 : ■ v 1 1 , header  ; 
if  lonil  then 
repeat 

1 inanum: *1 1 . 1 ine ; If  1 inenum>0  then  signold:*l  else  s1gnold:*-l; 

1 inenum: « s i gnol d*LTRAN[ A8S( 1 Inenum) ] ; 
wr ite( vertf  1  le, *  ' . 1 inenum;0) ; 

1  :  « 1 1 . nextl ine ; 
until  1  »ni 1 ; 

IF  verta  >  0  then  wr i teln( vertf lie) ; 
v  1  : *v 1 t . nextvert ; 
until  vl*ta11vert: 

{  do  Lines  } 

writeln(vertf11e,'L  *  \LCOUNT:0); 

1 1  :  -startl  ine; 
repeat 

L  INENUM  :»LTRAN[ll*.l  Inenum]  ; 

write(vertfi le , LINENUM : 0 , * : < ’ , VTRAN[1 1 t . Ini tvert] :0 . 1 . ' , VTRAN[ lit. final vart]:0. ’>: ' 
p : -1 1 t . header ; 
if  ponil  then 
repeat 

wri te( vertf lie, '( ' .pt.xp: 0:3, ', ' ,pt .yp: 0:3, ’)’); 
p :  *pt . nextpt; 
until  p  -  n  i  1  ; 
writeln(vertf11 3) ; 

1 1  :  *1 1  t  .  r.extl  Iro ; 
until  1 1  * taill Ine ; 

End;  {  WRITE  VERTICES  AND  EOGES  } 


{MAIN  PROGRAM  ST'.TTS 
Begin 

wr1teln( tty , 'Digltlxlng  Tablet;  Dutton  Functions:’); 

wr1teln(tty, ’  1:  start  a  new  vertex  (either  fresh,  or  as  the  end  of  the  currant'); 

writeln(tty. ’  arced  segment'); 

wr1teln( tty . '  2:  find  the  nsarost  vertex  to  this  point,  and  treat  It  as  the  end'); 

wr1t3lr.(tt,' .  ’  (or  otarr)  of  cjrror.t  srocd  5:c~:nt'); 


wr1teln(tty,  '  4:  taka  this  point'); 

wr1teln(tty , '  8:  lift  pen;'); 

wrltelnftty, '  If  nest  Is  1  or  2  then  CONTINUE  drawing'); 

wrltelnftty, ‘  If  nest  Is  4  than  do  a  DELETION'); 

wrltelnftty,'  If  nest  Is  8  then  QUIT’); 

DDGRAPH : -ASK ( Z  ) ; 

IF  ODGRAPH  then 
BEGIN 

DOGRAFX : -ASKf  3 ) ; 

GODGRAFX ; * ASKf  4) ; 

BOX FLAG: -ASKf  10) ; 

IF  80XFLAG  THEN 
Begin 

wrlteftty,'  Bos  s1ze(l .2) : ' ) ; 
readln(tty);read(TTY,boss1ze); 

End ; 

IF  GODGRAFX  THEN  GODINI: 

IF  DOGRAFX  THEN  ODINI; 

END; 

Rewr 1 te( vertf lie); 

f1rstp:*n11  ;  vptr : >ni 1 ;vertcount:-0;11necount:*0; 

INITLINELIST; INI TVERTLIST ; 

REPEAT  T:«  READINFILE  UNTIL  NOT ( T) ; 

TABLE TIO: -TRUE : DOTAB: -TRUE; 

IF  READIN  THEN  DOTAB : -ASK( 12 ) ; 

IF  DOTAB  THEN 
Begin 

wr 1 te( tty . ' TTY1 Ine  number:'); 
readln(  tty ) ; reao( tty , ttyl Ine) ; 
bi t Ini ( ttyl Ine) ; 

REPEAT  T : * INITTABLET  UNTIL  T;  {  get  corners  of  Image  } 

PPSET ; 

wrltelnftty. '‘Ready*' ); 

read1n:-false;f1nished:-false;toconf1rm:*fal se : jump : -true; LAST4 : -FALSE 

f Irs tp : *n1 1 ; 

repeat 

GETCOORDS; 
case  Bl  of 

'I' :STARTNEWVERTEX; 

' 2 ' : STARTOLDVERTEX ; 

’ 4 ' : NEWPOINT ; 

' 8 ' : Begin 

IF  NOT ( LAST4)  THEN 
Begin 

If  toconflrm  then  fin Ished : -true 
else  Begin 

toconf Irm: -true; 
wrltelnftty , '‘Pen  Is  UP*'); 

End ; 

Jump : -true 
End 

ELSE  wr1teln(tty. 'put  the  pen  down  first'); 

End; 

others -.Write  Inf  tty, 'Garbled  data’ ) 

end; 

until  finished; 

End; 

IF  OOGRAPH  THEN 
Begin 

IF  GODGRAFX  THEN  GOOF  INI  ELSE  IF  DDGRA'X  THE'I  DCFI.il; 

End ; 

WRITEFILE; 

End.  (MAIN) 
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:cod*  to  allow  PASCAL  to  gat  BITPAO  data  from  TTY  other  then  Job1*  TTY 
;  Macro-10  assembly  code,  with  local  code  for  TTY  driver 
title  trace 


1tem**0 

scl*-l 

SC2--2 

sc3**3 

sc4»»4 

p*-+017 

MODE--0  ; ASCII 
JOBFF*- tOt2 1 
chnstS*+07 16000000000 
ttyS(cp«rO04  7  0004001 16 
PP I 07 02 000000000 

DBEEP”  +  0047000400111  ;  BEEP  UUO  DEFN 

DPYCLR-' +0701 000000000 
INSKIP**+0051540000000 
SLEEP”+0047000000031 
OPDEF  PPSEL[PPIOE  0.] 

OPOEF  PPACTfPPIOT  1,] 

OPOEf  DPYPOS[PPIOT  2.] 

OPDEF  OPYS I Z[ PPIOT  3.] 

OPOEF  PPREL[PPIOT  4.] 

OPDEF  LEYPOS[PPIOT  6.] 

radl*  10 

entry  BITINI . B I  TP AD .PPSET.PPSET2, PPD0NE , TTY TAB 
BITINI:  POP  P.RETAO 

MOV EM  SC2 . LINE  ; LINE  PARAMETER  PASSED  BY  VALUE 

MOVE  I  SC3.+017 

BITIOP:  HRR 2  SC1.SC3  ; GET  A  CHAN  TO  TEST 

LSH  SCI, 18+5 

OR  SCI .[CHNSTS  SC2] 

XCT  SCI 

TRNE  SC2,  +0400000  -.USED  ALREADY? 

SOJGE  SC3 ,BI FLOP  ; TRY  UNTIL  NONE  THERE 

JIJMPL  SC3  ,  NOCHNS 

HRRZS  SC 3 

LSH  SC3.5 

HRLZM  SC3.CHAN  ;THIS  IS  THE  CHANNEL 

MOVE  I  TEM, [TTYSKP] 

I  OR  ITEM, CHAN 

MOVEM  ITEM, BITCH* 

MOVE  ITEM, LINE 

SETZ  SC2, 

LIKLOP:  IOIVI  ITEM. 10  ; TREAT  AS  DECIMAL  NUMBER  (ALTHOUGH  IT  IS  OCTAL) 

ADDI  SCI, +020 

LSHC  SCI, -6  ;  SIDE  IT  OVER 

JUMPN  ITEM, LINLOP  :GO  TILL  DONE 

HLRM  SC2 , 8 1 TDVC+1  : LINE  IN  SIXBIT 

BITOPN:  MOVE  ITEM, [OPEN  BITDVC] 

IOR  ITEM. CHAN  ;MAKE  OPEN  WORD 

XCT  ITEM 

HALT  BITOPN 

MOVE  I  SC1.  TTIBFS 

EXCH  SCI, JOBE* 

MOVE  I  TEM ,  r  I.  BUF  2] 

IOR  i TtM , CHIN 

XCT  I T  'M 

MOVE-*  'Cl.TTFF 


MOVE 

ITEM, [INPUT] 

IOR 

ITEM. CHAN 

MOVEM 

ITEM, INCMO 

JRST 

8RETAD 

NOCHNS: 

OUTSTR 

[ASCIZ/ 

No  channels  left.,  forget  It 

/] 

POPJ 

P, 

TTYSLP : 

MOVE  I 

0.0 

SLEEP 

0. 

TTYTAB: 

INSKIP 

0 

JRST 

DOBTPD 

SETOM 

1(P) 

POPJ 

P. 

DOB T PD. 

xcr 

BITCHK 

JRST 

TTYSLP 

PUSHJ 

P, BITPAD 

SETZM 

1(P) 

POPJ 

P. 

BITPAD: 

POP 

P.RETAD 

HRRZM 

SC2.RETVL1 

HRRZM 

SC3.RETVL2 

HRRZM 

SC4.RETVL3 

f  RC6: 

MOVE  I 

ITEM. 6 

MOVEM 

ITEM, COUNT 

PUSHJ 

P. GETBYT 

JRST 

FRC6 

MOVEM 

SCI. FSTVAL 

MOVE  I 

ITEM, 5 

MOVEM 

ITEM. COUNT 

PUSHJ 

P. GETBYT 

JRST 

FRC6 

INDON: 

MOVEM 

SC1.8RETVL3 

MOVE 

SCI, FSTVAL 

IOIVI 

SCI. +0100000 

MOVEM 

SCI  .8RETVL1 

MOVEM 

SC2.9RETVL2 

JRST 

8RETAD 

GETBYT : 

SETZ 

SCI. 

BITILP: 

SOSG 

IBF  +  2 

PUSHJ 

P.REFBUF 

I  LOB 

SC2.IBF+1 

CAIG 

SC2."9" 

CAIGE 

SC2.-0" 

JRST 

BITILP 

CNTLOP: 

IMULI 

SCI. 10 

ADDI 

SCl.-tO60(SC2) 

SOSG 

COUNT 

JRST 

COTCNT 

SOSG 

IBF+2 

PUSHJ 

P.REFBUF 

ILDB 

SC2.I8F+1 

CAIG 

SC2 , "9" 

CAIGE 

SC2 . "0" 

JRST 

3ITBAD 

JRST 

CNTLOP 

GOTCNT  : 

SCSJ 

IJF+2 

PUSHJ 

P.REFBUF 

ILDB 

SC2.I6F+1 

CAIG 

::2,"9" 

CAIGE 

JC2 . "0" 

JRST 

•- .OC.'.-D 

BAOCLR: 

SOSG 

t  F+2 

pi:  ■ ) 

\  :r.'T 

;  DONE  INITING  TTY  LINE 


; SLEEP  A  JIFFY 

;has  a  char  bean  typed? 

;  NO 

: return  TRUE  (SO  CALLER  GETS  INPUT  FROM  TTY) 


IF  THERES  A  CHAR  IN  BITPAD  BUFF.  DO  IT 
ELSE  SLEEP  A  BIT  AND  TRY  AGAIN 
GET  INPUT  FROM  BITPAD 
return  FALSE 


; FIRST  IS  6  DIGITS 


; INSIST  ON  GOOD  DATA 
: SECOND  IS  5  DIGITS 


; BAO  DATA.  LOOK  FOR  6  BEFOE  CONTINUING 
: RETURN  VALUE3  (Y  VALUE) 

; STRIP  OFF  BUTTON  FROM  FSTVAL 

; RETURN  BUTTON  (RETVL1) 

; RETURN  VALUE?  (X  VALUE) 


; RETURN  VALUE 
; COUNTER 


■OONT  START  TILL  GET  A  DIGIT 


; SUM  DIGITS 
; DIGIT  COUNTER 

;  BUFF  COUNTER 

; NOW  MUST  HAVE  5  OR  6  DIGITS  (1ST  OR  2ND) 


: GARBLED  ..  LOOK  FOR  6  AGAIN 
; GO  FOR  REST 
;3UFF  COL'. TER 

;nEXl  MUST  LE  NOI-DIGIT 


:  'ov  f:  j  cnr.ET 


— i  >n,  ,  t  .  Vi  -  -• 
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I  LOB 

SC2.I8F+1 

BITBAO:  CAIE 

SC2.Y016 

; CARRET 

JRST 

BAOCLR 

; NO ,  CLEAR  THE  CHARS 

MOVE  I 

SC1.-1 

OBEEP 

SCI. 

OUTSTR 

[ASCIZ/Bad 

data/] 

POPJ 

P. 

; YES .  RETURN  ERROR 

GOODE NO: AOS 

(P) 

-.GOOD  ONE  SKIPS  BACK 

POPJ 

P, 

REFBUF:  MOVE 

item.incmd 

XCT 

ITEM 

SOS 

IBF+2 

POPJ 

P, 

; GO  BACK  WITH  NEW  DATA 

PPSET :  PPSEL 

2 

; PAGE  PRINTER  2  ON 

PPACT 

tOlOOOOO 

;  SE  E  2,  eras*  0 

DPYPOS 

-450 

DPYSIZ 

f 0200 1 

;  2  LINES  AT  BOTTOM  OF  SCREEN.  SCROLL 

CONTINUOUSLY 

LEYPOS 

0 

;  LINE  EDITOR  AT  BOTTOM  OF  SCREEN 

PPACT 

t0500000 

.SEE  2  and  0  (now  cleared) 

POPJ 

P. 

PPSET2 :  PPSEL 

2 

;  PAGE  PRINTER  2  ON 

PPACT 

tOlOOOOO 

;  SEE  2.  erase  0 

DPYPOS 

-400 

DPYSIZ 

t0400 1 

;  2  LINES  AT  BOTTOM  OF  SCREEN.  SCROLL 

CONTINUOUSLY 

LEYPOS 

0 

; LINE  EDITOR  AT  BOTTOM  OF  SCREEN 

PPACT 

T0500000 

: SEE  2  and  0  (now  cleared) 

POPJ 

P. 

PPOONE.  PPACT 

t0400000 

; SEE  0  (clear  2) 

DPYCLR 

; FREE  PP  2 

POPJ 

P  . 

BITDVC :  EXP 

MODE 

SIXBIT/TTY  / 

EXP 

IBF 

IBF :  EXP 

0.0.0 

TTIBFS:  BLOCK 

100 

INTEGER 

INCMD, CHAN, LINE, RETAO.RETVL1.RETVL2.RETVL3. COUNT, FSTVAL.BITCHK 

END 

Epipolar  Program 


10:  Epipolar  Registration  Program 


See  section  4  for  documentation  on  using  this  program. 


Begin  "epipolar  determinations" 

REQUIRE  " 16A"  COMP ILER_SWITCHES ;  comment  FIXR.FLTR.ADJSP; 

DEFINE  w- ”3 . 141592663" : 

DEFINE  UNLESS*"WHILE 

DEFINE  ..-"COMMENT 

DEFINE  THRU*"STEP  1  UNTIL"; 

DEFINE  TAB*” ' 1 1 "  ; 

DEFINE  CR-"  ’  16"  ; 

DEFINE  IF-" ’ 12"  ; 

DEFINE  CRLF-"’16&'12"; 

DEFINE  CCRLF*"V15V12"; 

DEFINE  NN-20; 

REQUIRE  "CAMRAS . HDR[1 , jfm]"  SOURCE_FILE: 

a  Gennery's  camera  solver; 

REQUIRE  "MATRIX. HDRT1. JFM]"  SOURCE_FILE; 

n  Matrix  operation  utilities; 

REQUIRE  ”CHOSPT.HOR[S.HHB]"  SOURCE_FILE; 

o  Harlyn's  point  matcher; 

REQUIRE  "PIXSAIfS.HE]"  LIBRARY; 

REQUIRE  "PIXFAI[S.HE]"  LIBRARY: 

a  Source  files  for  image  manipulation  procedures:  see  below; 

EXTERNAL  INTEGER  PROCEDURE  P I XE L ( RE F E RE NCE  INTEGER  PIX;  INTEGER  I.J); 

o  gets  the  value  PIX  of  a  pixel  of  coordinates  I.J; 

EXTERNAL  PROCEDURE  PUTEL( RE F E RENCE  INTEGER  PIX;  INTEGER  1.0. VALUE); 

.>  sets  the  value  PIX  of  a  pixel  of  coordinates  I.J; 

EXTERNAL  INTEGER  PROCEDURE  PFLDIM( STRING  FILNAM); 

a  gets  the  proper  array  dimension  to  represent  the  image  called  FILNAM; 
EXTERNAL  INTEGER  PROCEDURE  GE TPF L ( STR ING  FILNAM;  REFERENCE  INTEGER  PICTURE): 
a  fills  array  starting  at  PICTURE  (  for  dimension  see  PFLD1M  ) 
with  picture  contained  in  file  FILNAM; 

EXTERNAL  INTEGER  PROCEDURE 

PUTPFL( REFERENCE  INTEGER  PICTURE;  STRING  FILNAM;  INTEGER  MODE(l)); 
q  writes  out  a  picture  file  FILNAM  with  data 
contained  in  array  starting  at  PICTURE; 

EXTERNAL  INTEGER  PROCEDURE  PI XDIM( INTEGER  HE IGHT , WIDTH , BITS) ; 

a  gets  the  proper  array  dimension  to  represent  a  HEIGHT  x  WIDTH  image; 
EXTERNAL  INTEGER  PROCEDURE  MAKP I X ( INTEGER  HE IGHT , WIDTH . BITS ; 

REFERENCE  INTEGER  PICTURE); 

a  creates  an  image  array,  starting  at  PICTURE,  for  a  HEIGHT  x  WIDTH  image; 
EXTERNAL  PROCECURE  WIPE( REFERENCE  INTEGER  PIX;  INTEGER  value): 
a  maxes  every  pixel  of  the  array  equal  to  value; 


a  Global  vari • bl _s ; 

REAL  ARRAY  GG[1:5]; 

a  Adjusted  values  of  the  parameters  ( azimuth . el evatlon , pan . ti 1 t , end  roll); 
REAL  ARRAY  KM[1:3.1:3]; 

a  rotation  matrix  between  basal  and  base  2; 

REAL  ALPHA. BETA. GAMMA , L AMOA . MU . NU ; 

a  translation  vector  components  in  base  1  and  2; 

INTEGER  ECASE; 

a  ECASE  is  a  flag  for  the  kind  of  geometry  we  have 
EC'SE *  1  :  two  epipoles  {  most  common  :"2  ) 
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ECASE-2  :  one  eploole  In  image  1 

ECASE-3  :  one  eplpole  in  image  2 

ECASE-4  :  no  epipole; 

INTEGER  PBITS; 

a  number  of  bits  per  pixel,  returned  by  CHOSPT: 

REAL  ARRAY  E 1 . E2[ 1 : 2 ] ; REAL  ARRAY  V1.V2[1:2]; 

a  points  and  vectors  in  images  1  &  2, 

epipoles  or  epipolar  directions  when  they  exist, 
else,  generic  points  and  vectors; 

REAL  F  F 1  ,  f  F  2  ; 

a  Focal  distances  in  images  1  &  2; 

REAL  ARRAY  M . MINV[ 1 : 2 . 1 : 2]  ; 

a  matrix  A  of  the  text  and  its  inverse; 

REAL  K.N1; 

a  aeneric  and  total  epiline  numbers:  OsK<Nl; 

REAL  L0X.L0Y.LNX.LNY; 

a  dimensions  of  the  images(  output  of  CHOSPT)  and  of  the  transforms; 
INTEGER  IK. ILOX . ILOY . ILNX . ILNY; 

o  the  same,  truncated,  as  integer  variables: 

REAL  THE TA. THE IA0 . THETA1  THE  TA2 ; 

a  gener i c . m i n i mum . max imum  and  increment 
angles  under  which  image  1  is  viewed  from  El; 

REAL  PH10.PHI1; 

a  min  and  max  angles  for  image  2; 

REAL  01MIN. DIMAX .Q2MAX .D2MIN.D; 

a  minimum,  maximum  and  generic  distances  from  epipoles  to  images  1  &  2; 
REAL  ARRAY  0F3[1 : 2 ] ; 

o  offset  for  case  3.  used  in  epiline  calculation; 

INTEGER  I;  u  generic  index; 

BOOLEAN  RIGHT1 . RIGHT2 ; 

o  true  if  epipole  right  of  image,  output  of  limits; 

BOOLEAN  OEBUG: 

STRING  P1.P2; 

a  filenames  for  images  PI  and  P2; 

INTEGER  SIZP.SIZS; 

a  Site  of  the  arrays  that  contain  original  and  transformed  image; 

REAL  ARRAY  L X [  1  : NN ]  ;  o  ARGUMENTS  FOR  CAMRAS; 

REAL  ARRAY  LY[ 1 : NN] : 

REAL  ARRAY  RX[ 1 : NM] ; 

REAL  ARRAY  RY[ 1 : NN] ; 

INTEGER  NUM; 


o  I/O  procedures; 

INTEGER  PROCEDURE  GET.V'SY'R; 

Begin  INTEGER  ANS;  AN5-INCHRW  LAND  '137;  PRINT (CRLF ) ;  RETURN(ANS)  End; 

BOOLEAN  PROCEDURE  ASK( STRING  QUESTION):  o  ask  a  yes/no  question; 

Begin  "ask" 

INTEGER  ASK; 

DO  BEGIN  pr int ( quest  ion ,"( Y  or  N)7  "  ) ; 
ask'-getinswer  ; 

END  UNTIL  ask  - "Y"  v  ask-MT; 
return(if  ask-”Y"  then  TRUE  else  FALSE) 

End  "  i  s  k "  ; 

PROCEDURE  GETD(REFERE  iCE  INTEGER  NUM;  STRING  QUESTION:  INTEGER  EFAULT) ; 
Begin  "GETO” 

PRINT] QUES  IC  ); 

PTOS  T  R  (  0  .  C\  -  .  ,  LLf )) ; 

NUM  -  CVC(  L) 


PROCEDURE  SPACEPAUSE( STRING  TEXT); 
BEGIN 

INTEGER  T; 

PRINT(TEXT); 

00  T-INCHRW  UNTIL  T-" 
PRINT(CRLE) ; 

END; 


PROCEDURE  GETO( REFERENCE  INTEGER  NUM;  STRING  QUESTION;  INTEGER  OEFAULT) ; 
Begin  "GETO" 

PRINT (QUESTION) ; 

PTOSTR ( 0 , CVOS( DEFAULT) ) ; 

NUM  -  CVO(  INCHWL ) 

End  "GETO"; 

PROCEDURE  GETF( REFERENCE  REAL  RNUM;  STRING  QUESTION;  REAL  DEFAULT); 

Begin  "GETF" 

STRING  NUMSIR; INTEGER  BRCHAR ; 

PRINT (QUESTION) ; 

PTOSTR(0 .CVF(OEFAULT) ) ; 

RNUM  -  REALSCAN(NUMSTR-INCHWL, BRCHAR); 

End  "GETF"; 


a  Arithmetic  procedures; 

SIMPLE  PROCEDURE  PR0D2  (REAL  ARRAY  A,8,C); 
a  multiplies  2*2  matrices; 

BEGIN 

REAL  T; 

INTEGER  I , J ; 

FOR  I-l  STEP  1  UNTIL  2  DO 
BEGIN 

T-0; 

FOR  J-l  ST E •  1  UNTIL  2  DO 
T  -  T  +  B[I,J]*C[J]; 

A[I]-T 

END 

END; 


SIMPLE  PRO.rDURE  PR0D3  (REAL  ARRAY  A,B,C); 
o  multiplies  3*3  matrices; 

BEGIN 

REAL  T; 
integer  i.j; 

FOR  1-1  STEP  1  UNTIL  3  DO 
BEGIN 
T-0; 

FOR  J-l  STEP  1  UNTIL  3  DO 
T  -  T  -  8[I . J]*C[J] ; 

A[I]-T 

ENO 

ENO; 


PROCEDURE  PARAMETERS; 

determines  the  rotation  matrix  R  and  the  translation  ’/actor  t  in  oth  bases; 
Begin  "PARAMETERS" 

REAL  ARRAY  TM1[ 1 : 3 . 1 ; 3] , T1 . T2[ 1 : 3] ; 

a  TH1  is  translation  ictri-.  (.-.rtrix  *  Ip  C  r;-) 
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T 1  and  T2  are  translation  vectors  In  base  1  and  base  2; 

SEAL  ARRAY  A1 . A2 , B 1 . B2 . B3 , D[ 1 : 3 . 1 : 3] ; 

Sl-S?-S3-3 ; 

R0TMAT(Al,0,GG[l].-2); 

R0TMAT(A2,O,GG[2],l); 

R0TMAT(B1,D.GG[3],2); 

fi0TMAT{B2.0,GG[4],-t); 

ROTMAT(B3,D,GG[5],-3); 

MULT ( TM1 , A1 , A2 ) ;  a  translation  matrix; 

MULT ( 0 , B2 , 81 ) : 

MULT ( RM , 83 , 0 ) ;  a  rotation  matrix; 

ALPHA*-T1[1]«-TM1[1,3];  a  translation  vector  is  3rd  column  of  translation  matrix 
BETA-T1[2]-TM1[2.3]; 

GAMMA*- T  1[3]*-TM1[3,3]; 

PR0D3 ( T 2 , RM , T 1 ) ;  a  T2  is  translation  vector  in  second  base; 

LAMDA*-T2  [  1  ]  ; 

MU*-T2[2]; 

NU-T2[3]; 

PRINT( "ROTATION  MATRIX  RM : " ,  CRLF ) ; 

PRINT(RM[1,1].RM[1.2].RM[1.3],CRLF); 

PRIN[(RM[2.1],RM[2.2],RM[2,3],CRLF); 

PRINT(RM[3, 1].RM[3,2].RM[3.3],CRLF); 

PRINT(" TRANSLATION  VECTOR  t:”.CRLF); 

PRINT(ALPHA. BETA. GAMMA. "IN  BASE  ONE”. CRLF); 

PRINT ( L AMDA , MU, NU . "IN  BASE  TWO", CRLF); 

end  "PARAMETERS"; 


a  Geometric  procedures; 

SIMPLE  PROCEDURE  RECT ( REAL  EX , EY , VX ,VY , RO ; REFERENCE  REAL  X,Y); 

a  rect  returns  the  rectangular  coordinates  of  a  point  with  olar  coordinates 
ro  and  a  with  respect  to  a  shifted  origin  e(ex,ey). 
vx  and  vy  are  supposed  proportional  to  cosa  and  sin  a,  but  vx>0: 
a  is  always  in  [-*/2,w/2],  and  RO  should  have  a  sign; 

BEGIN  "RECT" 

X-SQRT(VXt2+VYf2); 

VX*-VX/X ; 

VY*-VY/X: 

X*-EX*RO*VX ; 

Y*-EY+R0*VY; 

END  "RECT"; 


PROCEDURE  POLAR (  REAL  EX.EY.X.Y;  REFERENCE  REAL  THETA, RO); 

a  POLAR  returns  the  polar  coordinates  of  a  point  (x,y)  with  respect  to  a 
shifted  origin  (ex.ey) 

RO  is  always  a  positive  number, and  THETA  Is  In  the  determination  [0,2w]; 
Begin  "POLAR” 

RO-SQRTE(EX-X)t2+(EY-Y)t2); 

THETA*-ATAN((EY-Y)/(EX-X)); 

IF  ( X-EX ) <0  THEN  THETA-THETA+w ; 

IF  THETA  <0  THEN  THETA»THETA+2*w: 

END  "POLAR"; 
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SIMPLE  PROCEDURE  NEWX (  REAL  EX . EY . VX . VY ; REFERENCE  REAL  NX): 

a  A  point  E  has  the  coordinates  EX  and  EY  in  a  plane  where  there  1i  a  vector  V. 
If  this  plane  is  rotated  so  that  V  becomes  the  new  X-axis,  and  shifted  so 
that  the  old  origin  is  anywhere  on  the  new  Y-axis,  the  new  abscissa  of 
E  is  just  its  projection  onto  V:; 

BEGIN;  NX*-(  EX*  VX+E  Y*VY)  /SQRT(  VXT2+VY *2 ) ;  END ; 

BOOLEAN  PROCEDURE  INPICTURE ( REAL  X.Y.LX.LY); 

R£TURN(  IF  (QsXsLX  a  OsYsLY)  THEN  TRUE  ELSE  FALSE); 


PROCEDURE  L IMI TS(  REAL  LX  . LY , EX . EY ; REFERENCE  REAL  ALPHAO , ALPHA1 , OMIN.OMAX; 

REFERENCE  BOOLEAN  RIGHT); 

a  LIMITS  returns: 

1)  the  minimum  and  maximum  distances  OMIN  and  DMAX  from  a  point 

E  to  a  rectangle  of  length  LX  and  height  LY  with  bottom  left  corner  at 
origin. (  typically,  a  picture  and  its  epipole  ) 

2)  The  minimum  and  maximum  angle  through  which  the  rectangle  is  viewed 
from  point  E.  These  angle  ALPHAO  and  ALPHA1  are  in  the  determination 
[0.2w](  or  occasionally,  [-w/2 ,  -t-ir/2]  )  and  the  distances  are  consequently 
positive  only.  ; 

BEGIN  "LIMITS" 

REAL  X.Y; 

REAL  01, fi2.fi3.fi*'. 

DMIN*-SQRT (EXt2+EYt2);  a  initialize; 

DMAX**0 ; 

TOR  X«-l  THRU  LX  DO 
BEGIN 

«  DESCRIBE  HORIZONTAL  SIDES  OF  RECTANGLE; 

OMAX-OMAX  MAX  SQRT ((EX-X)t2+EY»2)  MAX  SQRT( ( EX-X ) * 2+( E Y-LY ) t2 ) ; 

DMIN-DMIN  MIN  SQRT(  ( EX-X  )  1 2*E Yt2  )  MIN  SQRT(  ( EX-X )  t 2  +  ( EY-LY )*2 ) ; 

END; 

FOR  Y-l  THRU  LY  DO 
BEGIN 

a  OESCRIBE  VERTICAL  SIDES  Of  RECTANGLE: 

DMAX*-DMAX  MAX  SQRT (  EX  t 2  +  (  EY-Y)  t2 )  MAX  SQRT(  ( EX-LX  )  tZ+(  E Y- Y) *2  )  ; 

DMIN-DMIN  MIN  SQRT( EXf2+( E Y- Y ) t2 )  MIN  SQRT( ( EX-LX ) r2*( EY-Y)*2 ) ; 

END; 

a  angle  determination; 

a  All  angles  are  initially  supposed  in  [0,2*]; 

IF (OsEXsLX  a  Os£YsLX)THEN  BEGIN 

a  epiline  is  inside  the  Image; 

ALPHAO-O; 

ALPHAl«-2*w; 

END 

ELSE  BEGIN 

a  el  1$  outside  the  Image; 

determine  the  limit  angles  by  Inspecting  the  corners  only; 

P0LAR(  EX  ,EY,0,0,/?1,D); 

POLAR(EX,EY,LX.O,02,D); 

POLAR(EX  ,EY,O,LY,03,O); 

POLAR (EX,EY,LX,LY,04,O): 

ALPHA1<-(01  MAX  02  MAX  03  MAX  04); 

ALPHAO-(01  MIN  02  MIN  03  MIN  04); 

IF  ALPHA 1-ALPHA0>»  THEN  BEGIN 

a  [ 0 , 2»]  Is  not  a  convenient  Interval  for  our  angles 

(  the  apipole  is  loft  of  the  Image,  on  the  same  vertical  level); 
a  Angles  are  In  [-»/2,*/2]; 

01*-ALPHA1 ; 

ALFHAli-ALDAO; 

ALPHAO*-01-2*w; 

IF  DEJUG  THEN  PR INT( "EPIPOLE  ON  LEFT  SIDE  AND  SAME  VERTICAL  LEVEL  AS  IMAGE ,CRLF) ; 
IF  DEBUG  THEN  PRI1)T("FCR  CONNEXITY,  ANGLES  ARE  NOT  CHOSEN  IN  [0.2*]"); 

I*  rEC'.'G  THEN  ?".::T("IN  LIMITS,  aO-  ”,  ALPHAO  ,CRLF )  ; 
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IF  DEBUG  THEN  PRINT("IN  LIMITS .  al-  " , ALPHA1 .CRLF ) ; 

END; 

a  if  epipole  right  of  picture,  sat  flag; 
RIGHT-(£X>LX/2); 

a  and  invert  scanning  order  by  exchanging  angles; 
IF  RIGHT  THEN  BEGIN 

01-ALPHA1; 

ALPHA1-ALPHA0; 

ALPHA0-/I1 ; 

END; 

END; 

END  "LIMITS"; 


PROCEDURE  FRAME (  REAL  EX , EY , VX . VY , LX , LV ; REFERENCE  REAL  I IX . I 1 Y . I2X , I2Y ; 
REFERENCE  BOOLEAN  OUTFRAME); 

a  Given  a  point  E  and  a  vector  V  (  they  define  a  line  ),  frame  will  determine, 
if  they  exist,  the  intersections  of  this  line  with  the  frame  of  the 
rectangle  [0:LX.O:LY]; 

BEGIN  "FRAME- 

REAL  X , Y ; INTEGER  FLAG; 


PROCEDURE  ASSIGN; 

BEGIN  "ASSIGN" 

IF  FLAG-0  THEN  BEGIN  I1X-X;  I1Y-Y;  END; 

IF  FLAG-1  THEN  BEGIN  I2X-X;  I2Y-Y;  END; 

IF  FLAG>1  THEN  PRINT( "ERROR  IN  PROCEDURE  FRAME" . CRLF ) ; 
FLAG-FLAG-1; 

END  "ASSIGN"; 


FLAG-0; 


IF  VX-0  THEN  BEGIN 
X-0  ; 

Y-EY+(X-EX)*VY/VX; 

IF  INPICTURE(X.Y.LX.LY)  THEN  ASSIGN; 
X-LX; 

Y-EY+(X-EX)*VY/VX: 

IF  INPICTURE(X, Y.LX.LY)  THEN  ASSIGN; 
END; 


IF  VY *0  THEN  BEGIN 
Y-0 : 

X-E  X- ( Y-EY) *VX/VY ; 

IF  INPICTURE(Y.X.LY.LX)  THEN  ASSIGN; 

Y-LY; 

X-E  X- ( Y-E Y) -VX/VY ; 

IF  INPICTURE(Y.X.LY.LX)  THEN  ASSIGN; 

END; 

OUTFRAME-(FLAG-O) ; 

IF  DEBUG  THEN  PRINT  ("  goes  through  :  ".I1X.I1Y,"  and  through:  " , I2X , I2Y.CRLF ) ; 
END  "FRAME"; 


PROCEDURE  OETCACE ; 

a  determines  the  nature  of  the  eplpolar  gjo,.,etry.  according  to  the  code: 
case  1:  two  epipoles  El  end  E2 
case  2:  one  epipole  El 
case  3:  one  epipolT  E2 
case  A;  r.D  .-pipolcd; 
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Begin  "OETCASE" 

REAL  THRESHOLD;  a  MINIMUM  POSITIVE  NUMBER; 

GETF(THRESHOLO, "MINIMUM  CONVERGENCE  ALLOWED:  ".0.0001); 

IF  ABS( GAMMA ) <THRESHOLO  THEN 

IF  ABS(NU)<THRESHOLO  THEN  ECASE*-4 
ELSE  ECASE-3 

ELSE  IF  ABS(NU)< THRESHOLD  THEN  ECASE-2 

ELSE  ECASE-1 ; 

PRINT("  CASE : " .ECASE .CRLF) ; 

and  "DETCAS £"; 


PROCEDURE  EPIPOLES; 

a  determines  the  eplpoles  and/or  the  aplpolar  directions 

It  also  calculates  some  parameters  used  by  EPILINE  (  matrix  M  and  offset  0F3  ); 


Begin  "EPIPOLES" 

REAL  OELTA;  a  Determinant  of  matrix  M; 

I F (ECASE - 1  v  ECASE-2)  THEN  BEGIN 

E1[1]-ALPHA*FF1/GAMMA; 

E1[2]-BETA*FF1/GAMMA; 

PRINT { "EPIPOLE  E1:",E1[1] ,  El[2] ,CRLF) ; 

ENO; 

IF(ECASE-1  V  ECASE-3)  THEN  BEGIN 

E2[l]-LAMOA*FF2/NU: 

E2[2]-MU‘FF2/NU: 

PRINT( "EPIPOLE  E2:".E2[1].E2[2],CRLF); 
o  M  Is  the  main  conversion  matrix; 
M[1.1]-NU*RM[1,1]-LAMDA*RM[3,1]; 

M[ 1 , 2]-NU*RM[ 1 . 2 ]-LAMDA*RM[3 . 2 ] ; 

M[  2 . 1 ]-NU*RM[ 2 . 1]-MU*RM[3 . 1] ; 

M[2,2]-NU*RM[2.2]-MU*RM[3,2]; 

a  Compute  the  inverse  of  matrix  M  for  transformation  P2-*P1; 
OELTA-M[ 1 , 1 ]*M[2 , 2]-M[l ,  2 ]*M[2 . 1 ] ; 

IF  OELTA-O  THEN  PRINT( "UNABLE  TO  INVERT  MATRIX  M  :  CAN'T  DUMP 

I")  ELSE 

BEGIN 

MINV[1,1]-  M[2.2]/DELTA; 

MINV[2 . 2  J-  M[1 , 1]/DELTA; 

MINV[1.2]-  -M[l , 2 ]/DELTA; 

MINV[2 ,1]-  -M[2 , 1 j/DELTA; 

END; 


IF  ECASE-3  THEN  BEGIN 

u  and  0F3  is  an  offset  used  In  ECArE  3; 
OF3[l]*-FFl*(NU»RM[l . 3]-LAMDA»RM[3 ,3]) ; 
OFSCE^FFl^CNU'RMfE.ShMU^RMCS.S]); 

END; 

END; 

IF(ECASE-3  v  ECASE -4)  THEN  BEGIN 

V1[1]-ALPHA; 

V1[2]-BETA; 

IF  V1[1]<0  THEN  BEGIN  Vl[l]— Vl[l];Vl[2;j— VJ[2]  Z  0; 
PRINT{ "EPIPOLAR  DIRECTION  VI : " . Vl[l] , Vl[2] .C.XF) ; 
END; 

IF(ECA3E-2  v  ECASE-4)  THEN  BEGIN 

V2[1]-LAMDA; 

V2[2]-MU; 

PR  I NT ("EPIPOLAR  DIRECTION  V2 : ” ,V2[1] ,VI[?],C.  IF); 
EKD; 
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End  "EPIPOLES"; 
a  more  procedures; 


PROCEDURE  EPIL1NE ; 

Begin  "EPILINE" ; 

IF  ECASE-1  THEN  BEGIN 

a  This  is  »n  epipolar  vector  in  PI; 

THETA«-THETAO+(  K+0 . 5)*THETA2 ;  a  take  middle  of  Interval; 

VI [  l]*-COS( THETA) ; 

V1[23*-SIN(  THETA) ; 
a  make  vl  point  to  the  right; 

IF  V1[1]<0  THEN  8EGIN  Vl[  1]— Vl[l  ] ;  V1[2]*--V1[2]  END; 
a  And  this  (V2)  is  the  corresponding  vector  in  P2; 

PROD2 ( V2 ,M . VI ) ; 

a  make  vl  and  v2  have  consistent  directions; 

IF  V1[1]*V2[1]<0  THEN  BEGIN  V2[l ]--V2[ 1] ; V2[2]— V2[2] ; END ; 


END; 


IF  ECASE-2  THEN  BEGIN 

a  This  is  an  epipolar  vector  in  PI; 

THETA-THETAO+( K+0 . 5 )*TH£TA2 ; 

Vl [  1 ]+COS( THETA) ; 

V1[2]+SIN(THETA) ; 

a  And  this  is  a  point  belonging  to  the  twin  line  in  P2; 
E2[1]-FF2*(RM[1,1]*V1[1]*RM[112]«V1[2])/(RH[3.1]»V1[1]+RM[3.2]»V1[2]); 
E2[2]-FF2«{RM[2.1]*Vl[l]+RM[2.2]‘Vlt2])/(RM[3,l]*Vl[l]+RM[3.2]»VI[2]); 
END; 

I F ( ECASE 1 3  v  ECASE-4)  THEN 

IF  ABS(Vl[2])*LOXsVl[l]*LOY  THEN  BEGIN 

a  We  pick  a  point  on  the  Y  axis; 

E  1[  1  ]+0 . 0 ; 

E1[2]-(IOX»ABS(V1[2]/V1[1])+LOY)»(K+0.6)/N1; 

IF  Vl[2]>-0  THEN  E 1[ 2 ]+E 1[2]-L0X»V1[2]/V1[1] ; 
ENO 

ELSE  BEGIN 

a  We  pick  a  point  on  the  X  axis; 
El[l]+(L0Y*ABS(Vl[l]/Vl[2])+L0X)»(K+0.6)/Nl; 

IF  Vl [2 ]>0  THEN  E  1[  1  ]+E  1[  1]-L0Y*V1[  1]/V1[2] ; 

E  l[2]+-0 .0; 

END; 

IF  ECASE-3  THEN  BEGIN 

a  V2  directs  the  twin  lijie  In  P2; 

PR0D2 ( V2 ,M, El ) ; 

V2[l]-V2[l]+OF3[l]; 

V2[2]«-V2[2]+0F3[2]; 

END; 

IF  ECASE-A  THEN  BEGIN 

a  The  twin  line  in  P2  goes  thru  E2; 
E2[1]+FF2»(RM[1.1]»E1[1]+RM[1,2]»E1[2]+RM[1,3]*FF1) 
/(RN[3.1]*E1[I]+R»»[3.2]»E1[2]+RM[3.3]*FF1); 
E2[2]«-FF2»(RM[2,1]*E1[1]+RH[2.2]*E1[*]+M[1.3]*FF1) 
/(RM[3,1]*E1[1]+RM[3.2]*E1[2]+RM[3,3]*FF1); 

END; 

E.td  -EPiL.' 
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Begin  "overht ad" ; 

a  0EBUG-ASK( "DEBUG");  DEBUG«-FALSE  ; 

IF  DEBUG  THEN  SETPR INT( "DEBUG" . "B”)  ; 

IF  DEBUG  THEN  SE  TFORMAT (0,7)  ELSE  SETFORMAT( 0 , 3 ) : 
CHOSPTS( LX , LY , RX . RY , NUM . ILOY . ILOX , PBI TS . PI , P2 , TRUE ) ; 
IF  DEBUG  THEN  PR  I NT ( "NUM-  -  ,  NUM ,CRLF ) ; 

IF  DEBUG  THEN  SETFORMAT( 0 . 7 )  ELSE  SETFORMAT(0 ,3) ; 
LOY*-ILOY;LOX«-ILOX; 

End  "overhead"; 


PROCEDURE  SOLVER; 

a  This  procedure  uses  Gennery's  code  to  do  the  camera  registration  from 
a  set  of  corresponding  points  in  the  two  images.  It  allows  the  user 
to  pick  the  points  on  the  images  in  an  interactive  way; 

Begin  "Solver" 


REAL  ARRAY  X 1 [ 1 ; num] . Yl[ 1 ; num] . X2[ I : num] , Y2[ 1 : num] ;  a  chosen  points; 
REAL  ARRAY  SXX[ 1 ; num] . SYY[ 1 : num] , $XY[ 1 : num] . 

GP[ 1:5], SP[ 1 : 5] ; 

REAL  ARRAY  S[l:5,l:6];  a  Covariance  matrix  of  the  errors  in  GG; 

REAL  ARRAY  RESID[ 1 : num. 1 ; 3] ; 

REAL  SD.ACC.Q.CONV; 

INTEGER  SDP.CORDST ; 

INTEGER  SDF.MAXEOIT; 

INTEGER  NREJECT ; 

BOOLEAN  ARRAY  RE JECT[ 1 : num] , FLAG[ 1 : 7] ; 

BOOLEAN  OUTPUT. RESOUT; 

FOR  1-1  THRU  NUM  DO 

Begin  "copy  values" 

Xl[I]«-LxjriJ;Yl[I)«-LYn]; 

X2[I]«-RX[I];Y2[I]*-RY[I]; 

End  "copy  values"; 

FOR  1-1  THRU  NUM  00 

SXX[I]-SXY[I].-SYY[I]-0; 

GETF(FF1,  "Focal  length  :",320);  FF2-FF1; 

GETO(MAXEDIT, "NumOer  of  points  that  can  be  re]9Cted : " , 0) ; 
SDP-0.67;CORDST-0;  SDF-0  ;  0UTPUT-RES0UT-( 1-1)  ; 
GP[l]«-1.67;GP[2]-GP[3]-GP[4]-GP[5]«-0.0; 


CAMERA( FF 1.FF2, SOP, CORDST.SDF, NUM. MAXEDIT. OUTPUT. RESOUT. LX. LY.RX.RY.SXX.SYY.SXY, 
GP, SP.GG.S.RESID. SO, ACC. Q.CONV, NREJECT. REJECT. FLAG) ;  o  do  the  solving; 

OUTSTR("Final  values  of  parameters  -”lcrlf); 

PRINT( "  azimuth:", CVF(GG[1])  ,crlf ) ; 

PRINT( "  elevation:" ,CVF(GGr 2] ).crlf); 

PRINT ( "  pan:",CVF(GG[3]).crlf); 

PRINT ( "  tilt:" ,CVF(GG[4] ) , crlf ) ; 

PRINT( "  roll :".CVF(GG[5]).crlf); 


End  "Solver"; 


PROCEDURE  NEWDIM; 

Begin  "New_dimensions" 

PRINT { CRLF ) ; 

GETF(N1. "NUMBER  OF  E?I?CL,.T  L;:;:S:  ".10); 
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CASE  ECASE  OF  BEGIN  "CASE  1" 

BEGIN  END; 

BEGIN 

LIMITS (LOX , LOY ,£2[1],E2[2],PHI0.PHI1 ,  02MIN , D2MAX .RIGHT? ) ; 

IF  INPICTURE  (£2[1].E2[2],  LOX  ,  LOY )  THEN  LNX-D2MAX+D2MIN 
ELSE  LNX-02MAX -02MIN ; 

LIMI TS( LOX , LOY ,E1[1].E1[2], THETAO , THETA1 , D1MIN , J1MAX , RIGHT1 ) ; 
a  theta?  is  the  incremental  angle; 

THE  TA2<-(  THETA1  -THETAO  )/Nl; 
a  the  epipole  is  either  inside  or  outside; 

IF  INPICTURE(E1[1].E1[-2].L0X.L0Y)  THEN  LNX-D1MAX+D1MIN  MAX  LNX 
ELSE  LNX-D1MAX-D1MIN  MAX  LNX; 

END; 


BEGIN 

NEWX ( LOX , LOY ,V2[1],V2[2],LNX) ; 

I  IMITS( LOX. LOY, E1[1].E1[2], THETAO. THETA1 , D1MIN . DIMAX . RIGHT1) ; 
IF  INPICTUR£(E1[1],E1[2],L0X,L0Y)  THEN  LNX-D1MAX+D1MIN  MAX  LNX 
ELSE  LNX-D1MAX-D1MIN  MAX  LNX; 

END; 


BEGIN 

NEWX(L0X.L0Y.V1[1].V1[2].LNX) ; 

LIMI TS( LOX. LOY, E2[ 1] ,  £2[2  ]  ,  PHI 0 ,  PHI  1 ,  D2MIN .  D2MAX .  RIGHT2 ) ; 

IF  INPICTURE(E2[l].E2[2].LOX,LOY)  THEN  LNX-D2MAX+D2MIN  MAX  LNX 
ELSE  LNX-DZMAX-D2MIN  MAX  LNX; 

END; 


BEGIN 

NEWX(L0X.L0Y,V2[1].V2[2].LNX)  ; 

NEWX( LOX . LOY ,V1[1],V1[2], LNY ) ; 

LNX-LNX  MAX  LNY; 

END; 

END  "CASE  1"; 

LNX-LNX+1  ; 

LNY-N1; 

a  prepare  the  integers  to  be  used  as  indexes  for  arrays; 

I  K-K  :  ILOX-LOX  ; ILOY-LOY  ;  ILNX»-LNX  ;  ILNY*-LNY; 

PRINT ( "THE  DIMENSION  OF  THE  NEW  PICTURES  WILL  BE:  ’ . ILNX . "( +4)  X  " . ILNY  ,  CRLF ) ; 
END  "New_dimensions" ; 


PROCEDURE  PRINTLINE; 

BEGIN  "PRINTLINE" 

BOOLEAN  OUTFRAME, MANUAL, PAUSE; 

REAL  X1.Y1.X2.Y2; 

INTEGER  IX1.IX2.IY1.IY2; 

print("  If  you  chose  the  step-by-step  option,  the  execution  stops  betwean" , crlf ) ; 
print(”  each  line,  waiting  for  you  to  type  Y",crlf); 

print("  If  you  do  not  chose  this  option,  the  execution  proceeds  automatically , " , crlf ) 
prir>t(”  either  in  SLOW  or  FAST  mode" . cr  1  f . crl f ) ; 

pr int( "When  printing  is  completed,  type  Y  in  both  cases  to  return  to  the  eiein  program 


MANUAL-  ASK( “  STEP-BY-STEP  PRINTING  "); 

IF  MANUAL  THEN  ELSE  PAUSE-ASX( "SLOW  MUTION"); 
INITOVERLAY ; 

FOR  K-0  THRU  Nl-T  DO 
BEGIN  "IVILOOP" 


IF  MANUAL  THEN  ASX("  LEFT  ") 

ELSE  IF  PAUSE  THEN  CALL  {1, "SLEEP"); 

EPILINE ; 

FRAME (El[l] ,El[2],Vl[l],Vl£2],LOX, LOY , XI , Y1 , X2 , Y2 .OUTFRAME) ; 
IF  OUTFRAME  THEN  PRINT( "EPILINE  #  ".K,"  OFF  PICTURE  l’.CRLF) 
ELSE  BEGIN 

IX1«-X1 ;  IX2--X2 ;  IY1«-Y1 ;  IY2«-Y2 ; 

SETLEFT; 

OVERLINE (IXt.IYl. 1X2, IY2 ) ; 

END; 

IF  MANUAL  THEN  ASX  ("RIGHT  ") 

ELSE  IF  PAUSE  THEN  CALL  (1, "SLEEP"); 
a  PRINT  (•'  In  picture  2,  "); 

FRAME ( E2[l) ,E2[2].V2[l],V2[2],LOX, LOY ,X1,Y1,X2,Y2 .OUTFRAME) ; 
IF  OUTFRAME  THEN  PRINT( "EPILINE  #  ”.X."  OFF  PICTURE  2",CRLF) 
ELSE  BEGIN 

SETRIGHT; 

IX1-X1 ;  IX2«-X2  ;  IY1«-Y1 ;  IY2«-Y2 ; 

OVERLINE ( 1X1, I Yl. 1X2, IY2); 

END; 

END  "MNLOOP" ; 

S PACE PAUSE ( "Type  <spece>  to  continue:  "); 

F INIOVERLAY ; 

END  "PRINTLINE"; 


PROCEDURE  NEWC1(  REAL  X , Y ; REFERENCE  REAL  Z ; REFERENCE  INTEGER  N); 

a  X  and  Y  are  the  coordinates  of  a  point  in  picture  2. 

NEWC2  returns  N  and  Z  which  are  the  number  of  the  epipolar  line(  new  ordinate) 
and  the  position  of  the  point  on  this  line  (  new  abscissa).; 

begin  "newel" 

IF  ECASE-1  THEN  BEGIN 

POLAR(  E1[1],E1[Z],X.Y, THETA, Z); 

I F  ( THETAOsTHETA+wSTHETAl )  THEN  THETA«-THETA«-w ; 

IF(  THETAOSTHETA-wSTHETAl )  THEN  THETA*-THETA-w; 

Z«-Z-D1MIN ; 

N«-N1*(  THETA-THETAO)  /  ( THETA1-THETA0)~0 . 5 : 

END  ELSE  PRINT( "SORRY,  DEGENERATE  CASE  NOT  YET  DEVELOPPED") ; 
end  "newel"; 

PROCEDURE  NEWC2(  REAL  X , Y ;REFERENCE  REAL  Z ; REFERENCE  INTEGER  N); 

a  NEWC2  does  the  same  as  newcZ  in  picture  *.  but  since  the  epilines  are  referenced 
by  their  position  in  picture  1.  this  Involve  some  transformation  first.; 

begin  "newc2" 

IF  ECASE-1  THEN  BEGIN 

POLAR ( E2[l] , E2[2] , X , Y , THETA ,Z) ; 

Z«-Z-D2MIN; 

V2[l]«-COS(THETA); 

V2[2]«-SIH( TnETA) ; 

a  And  this  (VI)  Is  the  corresponding  vector  In  PI; 

PROD2(Vl .MINV.V2) ; 

a  now  we  find  which  angle  v2  correspond  to  and  what  epiline  that  is; 
THETA-ATAN(V1[2]/V1[1]) ; 

IF( THET.*OsTHETA**ST!iETAl )  THEN  THETA*-TF:ETA«-w ; 
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I F  (  THETAOsTHETA-aSTHETAl )  THEN  THETA*-THETA-e ; 

N*-N1*(  THETA -THE  T  AO)/ (  THETAl-THETA0)-0.6; 

END  ELSE  PRINT ( "SORRY ,  DEGENERATE  CASE  NOT  YET  DEVELOPPED") ; 
and  "ne«rc2"; 


PROCEDURE  TRANSFORM; 

BEGIN  "TRANSFORM" 

BOOLEAN  OUTFRAME; 

REAL  X1.Y1.X2.Y2: 

INTEGER  IX1,IX2,IY1.IY2.NX,NN1 .VALUE ; 

SAFE  INTEGER  ARRAY  PP[0 ; SI2P-1 ] , SS[0 ; SI2S-1]  ;  a  INPUT  AND  OUTPUT  ARRAYS; 
INTEGER  XMIN.XMAX.N;  a  Length  info; 

STRING  EP1.EP2;  a  transforms  filenames; 

PROCEDURE  CLEARARRAY; 

BEGIN 

FOR  IK*-0  THRU  Nl-l  DO 
FOR  I-  0  THRU  ILNX-1  DO 
PUTEL ($S[ 0] . IK, I . 0) ; 

END; 


PROCEDURE  TRLINE (  INTEGER  IM;  REAL  EX , EY . DMAX . DMIN ; BOOLEAN  INVERT  ); 

BEGIN  "TRLINE" 

REAL  LOWERO . UPPERO ;  a  PICTURE  BOUNDS; 

REAL  VX.VY; 

IF  DEBUG  THEN  PRINT ( "EPILINE  »  " , K . CRLF )  ; 

EPILINE ; 

IF  IM~  1  THEN  BEGIN  VX-V1  [  1  ]  ;  VY«-V1[2]  ;  END 
ELSE  BEGIN  VX-V2 [  1  ]  ;  VY-V2 [ 2 ]  ;  END ; 

FRAME( EX. EY.VX.VY.LOX.LOY, XI, Y1.X2.Y2. OUTFRAME): 

IF  OUTFRAME  THEN  PRINT( "EPILINE  #  ",K."  OFF  PICTURE  ".IM.CRLF) 

ELSE  BEGIN 

I  K*-K  :  1 X  1«-X1 ;  I  X2*-X2  ;  IY1«-Y1 ;  IY2*Y2  ; 

IF  DEBUG  THEN  ELSE  OVERLINE ( IX I . IY1 . 1X2 . IY2 ) ; 

XMIN-DMAX; XMAX-0 ;  a  Initialize; 

IF  INVERT  THEN  BEGIN  UPPE RD— DMIN ; LOWERO— DMAX ; ENO 
ELSE  BEGIN  UPPE RD«-DMAX  ;  LOWERD«-DMIN ;  END; 

FOR  D<-  LOWERO  THHU  UPPERD  DO 
BEGIN  "pt-by-pt" 

RECT(EX.EY.VX.VY.D.Xl.Yl); 

IF  INP ICTURE ( X 1 , Y1 , LOX , LOY)  THEN 
BEGIN  "transfer" 

IXI*-X1 ;  IY1*Y1 ;  a  coordinates  on  source  image; 
VALUE«-PIXEL(PP[0],IY1.IX1)  ; 

IF  DEBUG  THEN  OVERLINE( 1X1 , IY1 , 1X1 , IYI ) ; 
NX-D-LOWERD; 

XMIN-XMIN  MIN  NX; 

XMAX-XMAX  MAX  NX; 

PUTEL (SS[0], IX. NX+4. VALUE); 

END  "transfer"; 

END  "pt-by-pt"; 


N*-XMI!I  LSH  -8;  a  8  MSBs  of  boginning  of  line; 

PUTEL (SS[0] . IX , 0 . N) ; 

N«-XMI,1  LAID  '377;  a  8  LSBs  of  beginning  of  line; 

PUTEL(SS[0],IK,1,N); 

N-XMAX  LSH  -8;  a  8  MSBs  of  end  of  line; 

PUTEL{ SSfOJ , IX . 2 , N) 

N*-Xf‘1X  LA  NO  '377; 

PUTEL'.cr^i.ir.a.”!) 


a  8  LSBs  of  end  of  line; 


a  Start  of  main  program; 

8£GIN"Main" 

OVERHEAD ; 

00  8EGIN 
SOLVER; 

PARAMETERS; 

DETCASE; 

EPIPOLES; 

NEWOIM; 

IF  ASK( "TRANSFORM  THE  IMAGES  INTO  EPIPOLAR  SPACE  ")  THEN 
IF  ECASE-I  THEN  BEGIN 

SIZP*-PFLDIM(P1); 

SIZP-PFLDIM{P2)  MAX  SIZP; 

ILNX«-I LNX+4 ;  a  make  space  for  line  length  info 
SIZS«-PIXDIM(  ILNY ,  ILNX  .PBITS) ; 

TRANSFORM; 

END 

ELSE  PRINT( "NO  TRANSFORMATION  FOR  DEGENERATE  CASE") 

ELSE  IF  ASK( "ORAN  EPIPOLAR  LINES  ")  THEN  PRINTLINE: 

IF  ASX("DUMP  THE  PARAMETERS  ")  THEN  DUMP; 

END  UNTIL  ASK( "ANOTHER  TRY  ") -FALSE ; 

END  "Main"; 


End  "epipolar  determinations"; 
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END: 


END  "TRLINE"; 


a  body  of  proc«dure  transform; 
a  D£BUG*-ASK( "DEBUG")  : 


PRINT("N£M  LEFT  PICTURE:  “): 

EPWNCHWL; 

PRINT( "NEW  RIGHT  PICTURE:  "); 

EP2«-INCHXIL; 

PRINT ( "PICTURE  I : " , CRLF ) : 

INITOVERLAY; 

SETLEFT; 

GETPFL(P1 ,PP[0)) ; 

MAKPIX(1LNY.1LNX.PBITS.SS[03)  ; 

FDR  K-0  THRU  Hl-1  00  TRL INE ( 1 , E If  1 ) . E IT 2 ] . DIMAX . D1MIN . RIGHT  1 ) ; 
PUTPFL(SS[0].EP1) ; 

FINIOVERLAY; 


INITOVERLAY; 

PRINT ( "PICTURE  2:". CRLF); 

SETRIGHT; 

GETPFL(P2,PP[0))  ; 

MAKPIX(ILNY.ILNX,PB1TS.SS[03): 

CUE ARARRAY ; 

FOR  K*-0  THRU  Nl-1  DO  TRL1NE(2  . E2[ 1)  . E2[2 3  .02MAX , D2MIN . RIGHT 2 ) ; 
PUTPFL(SS[0].EP2); 

SPACE  PAUSE (” Type  <space>  to  continue:  ”  ) : 

FINIOVERLAY. 

END  "TRANSFORM": 


PROCEDURE  DUMP; 
BEGIN  "DUMP" 


SE  TFORMAT (0.7); 

SETPRINT(NULL , "B" ) ; 

PRINT ( "ECASE *  ” , ECASE ,CRLF ) ; 

PRINT( "N1 *  ".Nl.CRLF); 

PRINT(-E1[1>  ".E1[1J.CRLF): 
PRINT(”E1[2]*  " . E 1[2] .CRLF  ) ; 
PRINT("E2[1]>  " . E2[ 1) , CRLF  ) ; 

PR INT ( "E2( 2 )»  " . E2[2 ] .CRLF  ) ; 
PRlNT("THETAO-  " , THE  TAO , CRLF  ) ; 

PR 1NT ( " THE  T A1 *  " , THETA1 .CRLF  ) ; 

PR INT ( "D1MIN*  * .01MIN.CRLF) ; 

PRINT ( "02MIN*  " , D2MIN.CRLF  ) ; 
PRINT("MINV[1.1)-  ".MINV[ 1,13, CRLF) 
PRINT ("MI NV[1,2]»  ",MINV[1,2],CRLF) 
PRlNT("MtNV[2.1]-  " . MINV[2 . 1] , CRLF ) 
PR1HT("MINV[2.23-  " ,MINV[2 . 2] .CRLF ) 

SETPR INT (NULL , "T" ) ; 

END  "CUVP" ; 
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